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ABSTRACT 

Aims. We compare the observed and theoretical parameters for the quiescent and outburst phases of the recurring nova T Pyx. 
Methods. IUE data were used to derive the disk luminosity and the mass accretion rate, and to exclude the presence of quasi-steady 
burning at the WD surface. XMM-NEWTON data were used to verify this conclusion. 

Results. By various methods, we obtained h djsk ~ 70 L G and M ~1.1 x 10~ 8 M yr~'. These values were about twice as high in 
the pre-1966-outburst epoch. This allowed the first direct estimate of the total mass accreted before outburst, M„ c „ =M pre -oB At, and 
its comparison with the critical ignition mass M, g „. We found M acc ,. and M jg „ to be in perfect agreement (with a value close to 5 x 
10~ 7 M o ) for M! ~ 1.37 M Q , which provides a confirmation of the thermonuclear runaway theory. The comparison of the observed 
parameters of the eruption phase, with the corresponding values in the grid of models by Yaron and collaborators, provides satisfactory 
agreement for values of Mi close to 1.35 M Q and logM between -8.0 and -7.0, but the observed value of the decay time 13 is higher 
than expected. The long duration of the optically thick phase during the recorded outbursts of T Pyx, a spectroscopic behavior typical 
of classical novae, and the persistence of P Cyg profiles, constrains the ejected mass M, g „ to within 10~ 5 - 10~ 4 M Q . Therefore, T Pyx 
ejects far more material than it has accreted, and the mass of the white dwarf will not increase to the Chandrasekhar limit as generally 
believed in recurrent novae. A detailed study based on the UV data excludes the possibility that T Pyx belongs to the class of the 
supersoft X-ray sources, as has been postulated. XMM-NEWTON observations have revealed a weak, hard source and confirmed this 
interpretation. 

Key words. Stars: novae - X-rays: binaries - Stars: supernovae: general 



1. Introduction 

Classical novae (CNe) are semi-detached binary systems in 
which a white dwarf (WD) primary star accretes hydrogen-rich 
matter, by means of an accretion disk, from a companion star 
which fills its Roche lobe. The "classical nova" outburst is a 
thermonuclear runaway (TNR) process on the surface of a white 
dwarf that is produced when, due to the gradual accumulation 
of hydrogen-rich material on its surface, the pressure at the bot- 
tom of the accreted layer, which is partially or fully degener- 
ate, becomes sufficiently high (P > 10 18 -10 2l) dyn cirT 2 ) for nu- 
clear ignition of hydrogen to begin. The critical mass for ignition 
Mig„ depends primarily on the mass of the white dwarf as M/ g „ 
~ M[ 7/3 , although more detailed models show that the ignition 
mass also depends on both the mass accretion rate M and the 
core WD temperature (Prialnik and Kovetz 1995; Townsley and 
Bildsten 2004; Yaron et al. 2005; see also Sect. 8). 

Recurrent novae (RNe) are a subclass of classical novae 
characterized by outbursts with recurrence time of the order of 
decades. We refer to Webbink et al. (hereinafter WLTO, 1987) as 
the seminal paper on the nature of recurrent novae and to Warner 
(1995), Hachisu and Kato (2001), and Shore (2008) for further 
considerations of this topic. Solid theoretical considerations in- 
dicate that a model of a recurrent outburst requires a high ac- 
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cretion rate M (10~ 8 -10 -7 M Q ) onto a WD of mass close to the 
Chandrasekhar mass limit (Starrfield 1985; WLTO 1987; Livio 
1994). However, nova models (Prialnik and Kovetz 1995; Yaron 
et al. 2005; see also Sect. 8) appear to evade these strict require- 
ments and allow the occurrence of recurrent outbursts also in a 
(slightly) less massive WD. The ejecta of RNe are expected to 
be less massive than those of CNe. This is because on the sur- 
face of the massive WD expected in a RN, the critical conditions 
for ignition are reached with a far less massive envelope. Studies 
of the ejecta of RNe indicate a mass between 10 6 -10 7 M Q , in- 
stead of the 10 4 -10 5 M observed in "classical" novae ( Gehrz 
et al. 1998, Hernanz and Jose 1998, Stan-field 1999). 

RNe represent a convenient laboratory to compare the pre- 
dictions of the TNR theory with the observations. From the ob- 
served M and the observed duration of the inter-outburst inter- 
val, in principle, one can obtain a direct estimate of the total 
mass accreted (M accr ) between two successive outbursts. This 
quantity can be compared with both the (theoretical) critical ig- 
nition mass Mi g „ and the mass of the ejected shell M c; , which 
can be estimated by spectroscopic methods (and, in principle, by 
photometric methods if dP/P can be accurately measured, Livio 
1991). A similar comparison cannot be made in the case of CNe 
because the amount of mass accreted prior to outburst is badly 
determined due to their long inter-outburst interval. 

Observations of RNe in quiescence and outburst can also 
be used to determine the secular balance between the total ac- 
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creted mass M accr and that ejected in the explosive phase M CJ , 
and therefore to investigate the possible role of RNe as progeni- 
tors of SN la. The ejecta of RNe have an almost "solar" chemical 
composition and are not enriched in "heavy" elements, such as 
carbon, oxygen and neon. This has been taken as an indication 
that the massive white dwarf in RNe is not eroded, will gain 
mass after each cycle of accretion and ejection, and will eventu- 
ally explode as a Type la supernova (Hachisu and Kato, 2002). 

In Gilmozzi and Selvelli (2007) (hereinafter Paper I), we 
studied the U V spectrum of T Pyx in detail and our main conclu- 
sion was that the spectral energy distribution (SED) is dominated 
by an accretion disk in the UV+opt+IR ranges, with a distribu- 
tion that, after correction for reddening E B „y=0.25, is described 
by a power law = 4.28 x 10~ 6 /T 233 erg cirT 2 s _1 A -1 , while 
the continuum in the UV range can also be well represented by a 
single blackbody of T ~ 34000 K. The observed UV continuum 
distribution of T Pyx has remained remarkably constant in both 
slope and intensity during 16 years of IUE observations. 

In the present study, we use this basic result from the IUE 
data and other considerations, to constrain the relevant parame- 
ters of the system and to attempt to understand the elusive nature 
of its outbursts. The five recorded outbursts of T Pyx occurred in 
1890, 1902, 1920, 1944, and 1966, with a mean recurrence time 
of 19+5.3 yrs (WLTO). All outbursts were similar in photomet- 
ric behavior and characterized by a decay time t3 ~ 90 d , a speed 
class that is substantially slower than in other RNe. We divide 
the paper into several self-contained sections. 

In Sect. 2, we discuss the problem of distance, which we 
determine to be 3500±500 pc. In Sect. 3, we analyze the mass 
function of the system. Using observational data and other con- 
siderations, we obtain 1.25 < M\ < 1.4 M Q , M 2 ~ 0.24 M Q , 
or M 2 ~ 0.12 M Q (depending on the as yet uncertain period), 
and i = 25 ± 5 degrees. In Sect. 4, we derive the absolute mag- 
nitude at minimum M v =2.53+ 0.23. In Sect. 5, we calculate the 
bolometric disk luminosity Ldi.sk ~ 70 L Q . 

In Sect. 6 we use various methods to determine the mass ac- 
cretion rate. The value depends on Mwd, and is M ~ 1.3 x 10~ 8 
M yr _1 for M W d = 1-33 M . In Sect. 7, we derive the accretion 
rate before the 1966 outburst, and find that it is about twice the 
value measured with IUE. In Sect. 8, we compare the accreted 
mass before the last outburst as determined from the M derived 
in Sect. 7, M accr ~ 4.5 xlO" 7 M Q , with theoretical calculations of 
the ignition mass. We find good agreement (including the inter- 
outburst interval) for M W d ~ 1-37 M G . In Sect. 9, we compare 
the observations of T Pyx with the theoretical models of Yaron 
et al (2005), finding good agreement for all parameters with the 
exception of f 3 . 

In Sect. 10, we reanalyze the historical data of the outbursts 
and reach the inescapable conclusion that the mass ejected in the 
explosion is M e j ~ 10~ 4 - 10 -5 M Q , about a factor 100 higher 
than the accreted mass. This discrepancy is discussed in detail 
in Sect. 11, and leads to the conclusion that the WD loses mass 
during the ourbursts. In Sect. 12, we revisit the parameters of the 
extended nebula around T Pyx, and demonstrate why the nebula 
cannot be used to derive information about the mass ejected dur- 
ing outburst. Section 13 discusses the mass balance and infers 
that T Pyx is not destined to become a SN la as often postu- 
lated for recurrent novae. Far from increasing its mass towards 
the Chandrasekhar limit, the WD in T Pyx is being eroded. 

In Sect. 14, we discuss the hypothesis that the WD in T Pyx 
is undergoing steady nuclear burning on its surface, and show 
that it is incompatible with both historical and UV observations. 
In Sect. 15, we report about XMM observations confirming our 
analysis by directly disproving the hypothesis. 



In Sect. 16, we discuss the recurrence time, and predict that 
the long awaited next outburst will occur around A.D. 2025. 
Section 17 reports our conclusions. 

2. The distance 

An accurate knowledge of distance is fundamental to the deter- 
mination of the disk luminosity ~Ldisk (and M„), the mass accre- 
tion rate M, and the mass of the observed nebula surrounding 
the system M„ e b (not to be confused with the mass ejected dur- 
ing outburst, M e j). 

In this section, we revisit the problem of determining the dis- 
tance to T Pyx. We adopt a new accurate determination of the 
reddening, E fi _y=0.25 ± 0.02, which is based on the depth of 
the A 2175 A interstellar absorption feature (see Paper I). 

An estimate of the absolute magnitude at maximum M™" 
can be obtained from the optical decay time t 2 or by extend- 
ing to a recurrent nova such as T Pyx, the applicability of the 
Maximum Magnitude versus Rate of Decline (MMRD) relations 
obtained using data of classical novae. 

This is justified by the fact that it is clearly established that 
the same kind of physical process for the outburst is involved i.e. 
a TNR on an accreting WD. The suggestion that recurrent novae 
in M3 1 are fainter than predicted by the MMRD curve (Delia 
Valle and Livio, 1998), applies only to the RNe subset that, un- 
like T Pyx, displays a short £3. The long tj, of T Pyx, the spectral 
behavior during outburst and the quantity of mass ejected (see 
Sect. 10) indicate that the same mechanism is operating as in 
CNe. 

The optical decay time of T Pyx is well determined (Mayall 
1967; Chincarini and Rosino 1969; Duerbeck 1987; WLTO 
1987; Warner 1995) with t 2 ~ 62 d and f 3 ~ 92 d . From these val- 
ues, we derived four different estimates for M"? ax (-7.01, -6.75, 
-7.00, -6.64) using the S-shaped MMRD relation of Delia Valle 
and Livio (1995), and the two linear and the one S-shaped rela- 
tions of Downes and Duerbeck (2000), respectively. 

An additional estimate of the distance can be obtained by 
assuming that T Pyx, as a moderately fast nova, radiated at near- 
Eddington luminosity at maximum. This provides 

LEddlLe = 4.6 x 10 4 (M, /M - 0.26) (1) 

For a massive white dwarf (Mi ~ 1.35 M ), this corresponds 
to M™* ~ -6.98. Close to the epoch of maximum light the ob- 
served color index showed oscillations about an average value of 
(B-V)=0.35 (Eggen et al. 1967, Eggen, 1968), or, after correc- 
tion for reddening, (B-V)„ ~ 0.10. Van den Bergh and Younger 
(1987) found that novae at maximum light have (B-V)„ = +0.23 
+ 0.06, while Warner (1976) reported that the intrinsic color of 
novae close to maximum is (B-V) fl ~ 0.0. The intermediate value 
for T Pyx was (B-V)„ ~ 0.10, which is indicative of a tempera- 
ture lower than 10,000 K and suggests that, close to maximum 
light, the nova radiated mostly in the optical. Since the bolomet- 
ric correction is quite small in this case (BC ~ -0.20), it follows 
that M" mx = -6.78, which is quite close to the values obtained 
above from the four MMRD relations. 

The mean value of the five estimates for M™ a t is -6.83, with 
a median value -6.79 and a formal standard deviation <x = 0. 16. 
The five estimates are intermediate between the value (-6.4) 
reported by Payne-Gaposchkin (1957), Catchpole (1968), and 
Warner(1995) and the single value (-7.47) given by Duerbeck 
(1981). 

Adopting M"? ax = -6.81, A v = 3.15 x E(B - V) = 0.79 
and m'l mx ~ 6.7 (Catchpole 1968, Warner 1995, Chincarini and 
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Rosino 1969), one obtains a distance d = lodHm-M-Ays) ~ 
3500 pc. The uncertainties in the individual parameters (0.16 for 
M™ ax , 0. 1 for m" MA , 0.06 for Ay) correspond to an uncertainty of 
about 330 pc in the distance, i.e. to a relative uncertainty of the 
order of 10 percent. 

A further estimate of the distance can be obtained by com- 
paring the dereddened flux of T Pyx at 1600 A with the accre- 
tion disk models by Wade and Hubeny (1998). As mentioned 
in Paper I, the most appropriate descriptions of the shape and 
depth of the lines (models ee and jj at low inclination angles) 
correspond to high values of white-dwarf mass and accretion 
rate, although they provide a continuum that is too steep. In 
any case, these high Mi, high M models, (we have extrapolated 
M to values as high as log M--7.5) clearly show that the ratio 
F mo l f F mo is about 100 °- Applying a scaling of a factor of ~ 
(1000) 5 = 31.5 to the model distance (100 pc) would recon- 
cile these two values. The derived value of ~ 3150 pc is in fair 
agreement with the assumed distance of 3500 pc. 

We recall that a range of distances to T Pyx can be found in 
the literature. Catchpole (1969) derived a lower limit of 1050 pc 
from the equivalent width of the Call K line (EW ~ 400 mA) and 
the EW-distance relationship derived by Beals and Oke (1953) 
for stars close to the Galactic plane. The EW measurement of 
Catchpole is however quite uncertain because of possible con- 
tamination by emission lines and the low quality of the spectra 
(the Call K line falls close to the edge of the spectrum and only 
two spectra were of sufficient signal in that region). 

Catchpole (1969), using the two spectra of the highest signal- 
to-noise-ratio found that the K line of calcium corresponds to 
radial velocity of +20 km s _1 . According to Catchpole (1969), 
this velocity agrees well with that of material travelling with the 
Galactic rotation in the direction of T Pyx; since this direction 
is, however, close to a node of the rotation curve, the Beals and 
Oke (1953) method cannot be used to determine the distance. 

Warner (1976), used the calibration given by Munch (1968) 
for stars close to a node in the Galactic rotation curve, and in- 
dicated that a star with a measurement of 0.4 A for the K line 
should be at least at a distance of 2 kpc. 

In a study of the IS Ca II K line towards O- and B-type stars 
in the Galactic disk, Hunter et al. (2006), presented a new cali- 
bration of the the total column density of Ca II K versus distance. 
In a plot of log N(K) versus log d, the estimated distance of 3500 
pc (log d = 3.545) would correspond to log N (Call) ~ 13.0. The 
EW of the K line measured by Catchpole (1969) corresponds to 
a column density of 4.9 X 10 12 (log N=12.69), which would cor- 
respond to a distance of 2200 pc. However, this value of column 
density is a lower limit, because the medium is not optically thin 
in the K line, and a yet higher value of distance is expected. 

By applying T Pyx theoretical predictions for light curves of 
CVs in outburst, Kato (1990) estimated a distance close to 4000 
pc, by comparing the observed magnitude at outburst with the 
theoretical absolute magnitude prediction at outburst. 

In the following, we assume a distance of d = 3500 + 350 
pc. A discussion of the critical role played by distance in the 
interpretation of the nature of T Pyx is given in Sect. 11.1. 



3. The system parameters 

We aim to determine both the disk luminosity and the mass ac- 
cretion rate of T Pyx. This requires prior knowledge of the sys- 
tem inclination angle i and the mass of the primary Mi. In the 
disk geometry, the inclination angle is critical to the estimate 
of the disk luminosity, while Mi and Ri (a function of Mi) are 



key parameters in the correlation between M and L^. For most 
CVs, the determination of these and other parameters (M2, P or b, 
2K) entering the mass function: 



(M 2 ■ sinif 
(Mi + M 2 ) 2 



= 1.037 x 10~ 7 -K\ P 



(2) 



is a difficult task, and accurate solutions have been obtained only 
for a few eclipsing systems. 

For T Pyx, the situation, "prima facie", does not appear en- 
couraging because no system parameter (not even the orbital pe- 
riod) have been accurately measured. Therefore, one has to start 
from a restricted range of values corresponding to the most accu- 
rately known parameters in the mass function equation, based on 
theoretical or semi-empirical considerations, and derive the cor- 
responding range of allowed solutions for the unknown relevant 
parameters. 



3. 1 . The orbital pehod(s) and K\ 

There have been several attempts to measure the orbital period 
of T Pyx, although, disappointingly, they have resulted in a wide 
range of values (see Paper I); this is a clear indication that there 
is no definite photometric clock associated with the system. 

Patterson et al. (1998) detected a stable photometric wave at 
P A =l /, .829 ± A .002 (a value close to that of the "most proba- 
ble" photometric period P h = 1 ;, .828 of Schaefer et al., 1992), 
although this interpretation is inconsistent with the presence of 
another signal at 2 A .635. 

The only spectroscopic determination of the period is that 
by Vogt et al. (1990), who, in a preliminary study based on a 
large number (101) of spectra of limited resolution (~3 A), re- 
ported a spectroscopic period of P/, ~ 3\439, without quoting 
uncertainties. This spectroscopic period is almost twice the pho- 
tometric value determined by Szkody and Feinswog (1988) in 
the infrared. In other CVs, this behavior is generally interpreted 
as evidence of ellipsoidal variations in the companion, in a sys- 
tem observed at a high inclination angle. This is not the case for 
T Pyx because, as convincingly shown in Paper I, the companion 
does not contribute to the J light curve, unless the distance is less 
than 200 pc. 

In the following, we consider system solutions for both P/, 
= l\829 (photometric, Patterson et al, 1998) and P/,=3\439 
(spectroscopic, Vogt et al. 1990). The inclusion of the spectro- 
scopic period could be criticized since it is the result of a prelimi- 
nary analysis that has not been confirmed by subsequent studies. 
However, given the complex photometric period structure and 
in the absence of a definite physical interpretation of the stable 
signal with P/ ! =l' ! .829 (that could even arise from the rotation 
of a magnetic white dwarf), we believe that the spectroscopic re- 
sults should be considered, since, in any case, they have not been 
challenged by similar studies. 

Vogt et al. (1990) found that the projected orbital velocity 
Ki of the primary, derived from an analysis of radial velocities 
of the emission lines, is approximately 24+5 kms -1 . Taking into 
account the limited resolution and the large number of spectra, 
this value is probably an upper limit. We assume that this value 
of Ki is representative of the primary star radial velocity (in- 
dependently of the period being considered), although we are 
aware of the problems encountered in associating radial velocity 
changes with the motion of the WD. 
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3.2. The mass of the primary 

Both theoretical considerations and observational determina- 
tions (although limited) indicate that the WD masses in classical 
nova systems are about 1.0 M Q , i.e. higher than inferred from 
the standard field WD mass distribution (Ritter et al., 1991). 
Theoretical models of recurrent novae require an even more mas- 
sive white dwarf (WLTO, Shara 1989, Livio 1994). Therefore, in 
the following, we allow the mass of the primary, Mi, to vary be- 
tween 1.25 and 1.4 M . 

3.3. The mass of the secondary 

In the absence of any direct information about the secondary 
(the accretion disk is the dominant luminosity source from the 
UV to the infrared, Paper I), we can estimate its mass M2 on the 
basis of the commonly accepted assumptions that: a) the sec- 
ondary fills its Roche lobe, b) its radius R2 is determined by the 
Roche geometry, and c) the secondary star obeys a mass-radius 
relation valid for low mass main sequence stars. In this standard 
description, the radius of the secondary is well approximated by 
the equivalent radius R^ of the Roche lobe itself, i.e. R2 ~ Rz. 
and the spherical volume of the secondary is assumed to equal 
the volume of its Roche lobe. The equivalent volume radii of the 
Roche lobe are given in tabular form by several authors (Plavec 
and Kratochvil 1964; Kopal 1972; Eggleton 1983; Mochnacky 
1984). For our purposes, it is convenient however to use the an- 
alytical approximation by Paczinsky (1971), which is valid for 
M 2 /Mi < 0.8: 

R L /A = 0.457256(M 2 /(Mi + M 2 )) I/3 (3) 

where we have introduced a new value for the constant (in- 
stead of 0.46224) because it yields results that differ by less than 
1% from those tabulated by Eggleton (1983 ) and Mochnacky 
(1984). 

By combining equation[3]with Kepler's third law: 

A = 0.50557 'M^Pjf 3 (4) 

one obtains an approximate relation between P, M2, and R2 (~ 
Rl): 

P h = 8.997/?2 /2 M 2 1/2 (5) 

In the crude assumption M2 ~ R2, this provides an approximate 
relation: M2 = 0.1 1 1 P/, by which a rough estimate of the mass 
of the secondary can be obtained if the orbital period is known. 
In general, however, the R-M relation for the low MS stars is 
described more accurately by a relation of the form: 

R 2 ~/3M% (6) 

By combining eq.|5]and eq[6j we derive a more general relation 
between M2 and P/, : 

Ml _ p(2/(3a-l)) _ g 997 (2/(l-3 ff )) _^(3/(l-3a)) 

Values for the parameters a and ft of the low mass main- 
sequence stars were derived from observational data and theo- 
retical considerations. In particular, Warner (1995) assumed that 
R2=M 2 3 ^ 15 (which provides an approximate fit to the data set by 
Webbink, 1990) and derived M^O.OeSxP 1 25 , while Smith and 
Dhillon (1998) found a constitutive relation R 2 =0.91 M° 75 and 
derived M=0.038xP 158 . 



Table 1. M 2 for P^l.829 and PM.439 from various M 2 -P re- 
lations or different a and /3 values in the literature (see text). 





P,,= 1.829 


P,,=3.439 


Patterson et al. 


0.134 


0.234 


Smith and Dhillon 


0.100 


0.267 


Knigge 


0.140 


0.236 


Warner 


0.138 


0.304 


Equation (7) 


0.121 


0.227 



Patterson et al. (2005) derived an empirical mass-radius se- 
quence for CV secondaries based on masses and radii mea- 
sured primarily for the superhumping CVs. They found that 
R 2 =0.62xM^ 61 for P shorter than 2.3 hours, and R 2 =0.92xM° 71 
for P longer than 2.5 hours, and derived M 2 =0.032xP 238 (p' < 
2.3) and M 2 =0.026xP L78 (P > 2.5). We inserted the a and (3 val- 
ues of Patterson et al.(2005) in Eq.|7]and obtained slightly differ- 
ent results, that is, M 2 =0.02823 P 2 41 (short P) and M 2 =0.02554 
P 1 77 (long P). 

Knigge (2006) obtained an independent M2-R2 relation by 
revisiting and updating the mass-radius relationship for CV sec- 
ondaries determined by Patterson et al. (2005). We have fitted 
with a power-law the data contained in his Table 3 and derived 
M2-P relations for short P and long P systems. The relations are 

M 2 = 0.0371 x P 2 - 2022 (P < 2.2), (8) 
M 2 = 0.02024 x P 1 ' 9802 (P > 3.2). 

Table 1 provides the M2 values obtained by these various 
methods, for P/,=l' ! .829 (photometric) and P/,=3''.439 (spectro- 
scopic). 

In the following, we therefore study the system solutions for 
M 2 = 0.12 +0.02 M and M 2 = 0.24 + 0.02 M , separately. 
The corresponding spectral types are in the range M3V-M5V 
(Kirkpatrick et al. 1991, Legget 1992). 

3.4. The system inclination 

The above derived range of allowed values for the system pa- 
rameters (Mi: 1.25-1.4 M Q , M 2 = 0.12 M Q (for P,,=l /! .829), M 2 
= 0.24 M (for P/,= 3'\439), Ki: 24 +5 km s" 1 ) can be inserted 
into Eq. (2) to find the corresponding range of solutions for the 
system inclination, which is considered to be a free parameter. 

Table 2 clearly shows that, given the observational and the- 
oretical constrains for P, Ki, Mi, and M2, the solutions for the 
inclination i are in the range between 20 and 30 degrees, a value 
close to 25 degrees being the most likely value. In particular, 
for Mi ~ 1 .35 M , (see also the following) the solutions for i 
do not change significantly, and for P /i =l /, .829 (M 2 =0.12 M Q ) 
and P A =3\439 (M 2 =0.24 M ) are close to 30 and 20 degrees, 
respectively (see also Fig. 1). 

Patterson et al. (1998) suggested that T Pyx is being observed 
at low system inclination, i.e. i ~10°-20°, due to the low ampli- 
tude of the orbital signal. Shabbaz et al. (1997) estimated the 
binary inclination i of T Pyx by measuring the separation of the 
H a emission-line peaks, and obtained a lower limit of ~ 6°. A 
low system inclination is also consistent with the small value 
of the radial velocity, the sharpness of the emission lines in the 
optical (Warner, 1995), and the steepness of the UV and optical 
continua (Paper I). However, the presence of radial velocity vari- 
ations (Vogt et al. 1990), and the modulations in the photometric 
variations preclude an i value close to zero. 
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Table 2. The system inclination for M, = 1.25-1.40 M and Ki 
=24 + 5 km s- 1 for the cases (1) P /l =l /l .829, M 2 =0.12; and (2) 
? h =3 h A39, M 2 =0.24 



Mi 


K, 


( (P=l".83) 


; - (P=3".44) 


(M Q ) 


(km s^ 1 ) 


(degrees) 


1.25 


19 


22.9 


14.7 


1.25 


24 


29.4 


18.7 


1.25 


29 


36.4 


22.8 


1.30 


19 


23.5 


15.0 


1.30 


24 


30.2 


19.1 


1.30 


29 


37.4 


23.3 


1.35 


19 


24.0 


15.4 


1.35 


24 


31.0 


19.6 


1.35 


29 


38.4 


23.9 


1.40 


19 


24.6 


15.7 


1.40 


24 


31.7 


20.0 


1.40 


29 


39.5 


24.4 




0.8 1 1.2 1.4 1.6 1.8 1 1.2 1.4 1.6 1.8 

M, [MJ M, [M ] 



Fig.l. The q vs Mi plane for P/,=3\439 and P /l =l''.829. The 
lines of constant inclination refer to K\ = 19 km s _1 (solid) and 
K\ = 29 km s _1 (dashed), spanning the error range in the value 
of K\ . The values of ;' are indicated. The error bars represent the 
ranges in M\ and M 2 discussed in the text. 



It should therefore be noted that, in spite of the large un- 
certainties in the observational data, the adoption of plausible 
theoretical assumptions and some semi-empirical constraints has 
enabled a quite restricted range for the values of the system incli- 
nation to be defined. In the following we will assume i = 25 ± 5 
degrees for the system inclination angle. 



4. The absolute magnitude at minimum 

During the quiescent phase after the 1967 outburst, the aver- 
age optical magnitude of T Pyx is m v , = 15.3 + 0.05 (Warner 
1995, Duerbeck 1987, Schaefer 1992); for d = 3500 + 350 pc 
and Ay=0.79 + 0.06, we obtain an (observed) absolute magni- 
tude at minimum M° b ' 5 =1.79 + 0.21 (Warner 1995 defined this 
to be the apparent absolute magnitude since the observed flux in 
quiescence depends on the inclination angle). 

Assuming that the visual radiation originates in a non- 
irradiated disk (Paper I), we correct the apparent absolute mag- 
nitude M° hs (i) for the inclination angle of the disk (see Warner 



1987, 1995) by a term AM v (z) to obtain the "standard" or "refer- 
ence" absolute magnitude: 

M corr = M ob S{i) _ AMi (/) (9) 

where AM v (z) = -2.5 log(2 cos ;) according to WLTO (1987), or 
AM v (i) = -2.51og[(l + 1 .5 cos cos i] according to Paczynski 
and Schvarzenberg-Czerny (1980). 

For ; ~ 25°, one obtains a correction of AM v (i) ~ 0.74 
± 0.06, by averaging these two relations, and M c ° rr = +2.53 + 
0.22 for the absolute magnitude, averaged over all aspect angles. 
This value is about 1.4 magnitudes brighter than the mean ab- 
solute magnitude of nova remnants in the same speed class (see 
Fig. 2.20 and Table 4.6 of Warner 1995). It is also about 2.0 
magnitudes brighter than the mean absolute magnitude of no- 
vae at minimum, as obtained from the values given in Table 6 of 
Downes and Duerbeck (2000). 

The T Pyx remnant is therefore one of the brightest nova 
remnants, its absolute magnitude My°"=+2.53 being close to 
that of the ex-nova HR Del, whose corrected value, M™"=2.30, 
implies that it is a quite luminous object (see also Paper I). 

5. The disk luminosity from IUE and optical data 

As shown in Paper I, the integrated UV continuum flux of T 
Pyx in the wavelength range 1 180-3230 A , after correction for 
reddening, is 1.94 x 10~ 10 erg crrrV 1 . For d = 3500 + 350 pc, 
the corresponding integrated luminosity of the UV continuum is: 

L uv ~ 2.85 x 10 35 erg s' 1 ~ 74.2 ± 15.0 L Q . 

where the 20 percent is caused by the 10 percent uncertainty in 
the distance. Most old novae have Lyy in the range 1-20 L Q , 
with V 841 Oph at ~ 24 L G (Gilmozzi et al., 1994), and V603 
Aql and RR Pic at ~ 1 1 L Q and ~ 3 L Q , respectively (Selvelli, 
2004). Therefore, T Pyx is also quite bright in the UV and, again, 
is challenged only by HR Del that has Luv = 56 L G (Selvelli and 
Friedjung (2002). 

In the following, we assume that the observed UV and op- 
tical luminosity arise from an accretion disk heated by viscous 
dissipation of gravitational energy. This agrees with the behav- 
ior of similar objects, like old novae and nova-like stars. This 
assumption is also strongly supported by the fact that the con- 
tinuum energy distribution in the UV range is well reproduced 
(Paper I) by a power-law with spectral index A ~ 2 33 , as predicted 
by theoretical optically thick accretion disk models radiating as 
a sum of blackbodies. 

The bolometric disk luminosity Ldisk can be estimated from 
the observed UV and optical luminosity, where the bulk of the 
continuum radiation is emitted, after correction for the inclina- 
tion and the unseen luminosity in both the infrared and at A < 
1200 A. The radiation emitted at wavelengths shorter than Ly ff is 
strongly absorbed and the energy is redistributed to longer wave- 
lengths (Nofar, Shaviv and Wehrse 1992, Wade and Hubeny 
1998). 

The integration (from A ~ 1000 A to the IR) of the power- 
law distribution which represents the observed UV and optical 
continuum of T Pyx, corresponds to a flux of about 3.6 x 10~ 10 
erg cirT 2 s _1 , which corresponds to a bolometric luminosity of 
5.24 X 10 35 erg s _1 ~ 136.5 L . This value refers to the adopted 
inclination of about 25 degrees. After correction to the "stan- 
dard" inclination of about 57 degrees and by considering an aver- 
age limb-darkening factor similar to that in the optical, we obtain 
an angle (47r) averaged bolometric disk luminosity (hereinafter 
Lrfijjt) of about 70+15 L , where the relative uncertainty (21 
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percent) derives from the combination of the uncertainties in the 
distance and the inclination. We consider this to be the reference 
disk luminosity for the "recent", post-1967-outburst epoch. We 
note that the Wade and Hubeny (1998) models indicate (their 
Fig. 8) that in the UV region close to 1200 A (where the effect 
is larger) the limb-darkening correction, from i of about 25 de- 
grees to i of about 60 degrees, corresponds to a flux reduction of 
about 30 percent (this is limb darkening alone, and the geometric 
projection factor cos ; is not included). 

A rough estimate of the An averaged bolometric disk lumi- 
nosity can also be obtained from M € l orr = +2.53. Adopting an 
average disk temperature of about 28,000 K (the UV continuum 
indicates ~ 34,000 K, although this is too hot to match the optical 
fluxes), the corresponding bolometric correction is B.C. ~ -2.7, 
and Lboi ~ 92.9 L e , in fair agreement with the previous estimate. 
In this comparison, we neglected the effects of temperature strat- 
ification in the accretion disk. 

6. The mass accretion rate 

One could in principle estimate the mass accretion rate M by 
comparing the observed spectral distribution with that of proper 
accretion disk models. This approach is not viable in our case 
for the following reasons: a) the spectral index of model disks 
that consist of blackbodies is hardly sensitive to the mass accre- 
tion rate M; b) most models providing M as a function of the 
continuum slope do not include the case of a massive WD with 
high accretion rates; c) disk models depend on a large number 
of parameters, so that model fitting does not generally provide 
unique results. 

Alternatively, if, as in our case, one can make the reasonable 
assumption that the disk is heated by viscous dissipation of grav- 
itational energy, the mass accretion rate M can be obtained from 
the relation: 

M = (2R\L disk )/(GM\) (10) 

With this method, M is not model dependent but the knowl- 
edge of Ldisi and M x is required. 
Numerically, M can be represented by: 

M = 5.23 x 10" 1 °(f>L disk /L Q with (11) 

= (R x /M x )/(R Xo /M Xo ) = 0.1235 • R x /M x , (12) 

where Ri is the WD radius in 10 3 R Q , Ri c = 8.10 x 10 3 R Q 
is the radius of a WD of mass Mi = 1.0 M Q , and Mi is in solar 
masses. 

We obtained average values for Ri as a function of Mi from 
various WD radius-mass relations in the literature (Hamada and 
Salpeter 1961; Nauenberg 1972; Anderson 1988; Politano et al. 
1990; Livio 1994). By fitting a quadratic function to these aver- 
age values, for Mi in the range 1.0 to 1.4 M Q , we found that: 

R X /R Q = -0.01315M 2 + 0.01777M, + 0.00347. (13) 
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The upper curve in Fig. 2 is a plot of the standard Ri ~ M l 
relation, which is generally valid for Mi < 1.0 M , while the 
lower curve represents the quadratic fit to the average values for 
Mi in the range 1.0 to 1.4 M Q , with a relative uncertainty of 
a Rl /R x = 0.07. 

Columns 2, 3, and 4 in Table 3 indicate Ri, <p, and M for M x in 
the range 1.00 - 1.4 M G . We note that <p = 1.00 for Mi = 1.0 M Q , 
and decreases to about 0.23 for Mi =1.4 M Q . The accretion rate 
M in Table 3 is calculated for the (present) luminosity ~ 



T 




2 - 



1 1.1 1.2 1.3 1.4 

Mass M, [MJ 

Fig. 2. The lower curve represents a quadratic fit to the average 
values of the WD radius as a function of Mi derived from vari- 
ous WD radius-mass relations in the literature (dots). The upper 
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curve is the standard Ri ~ Mj relation. 

Table 3. Mi, Ri, the factor cp =(Ri/Mi)/(Ri /Mi ), the post- 
1967 mass accretion rate M (for l^disk = 70 L ), the ignition 
mass, and the recurrence time T=M ign /M (see Sect. 8). 



Mi 






M 




T 


(Mo) 


(10- 3 R o ) 




(10- 8 M Q yr 


') (10- 7 M o ) 


(yrs) 


1.00 


8.10 


1.000 


3.66 


108.78 


297.2 


1.05 


7.63 


0.897 


3.28 


92.88 


283.2 


1.10 


7.10 


0.797 


2.92 


77.75 


266.3 


1.15 


6.51 


0.699 


2.56 


63.32 


247.3 


1.20 


5.85 


0.602 


2.20 


49.19 


223.6 


1.25 


5.14 


0.508 


1.86 


35.21 


189.3 


1.30 


4.35 


0.413 


1.51 


22.22 


147.1 


1.33 


3.85 


0.357 


1.31 


14.78 


112.8 


1.35 


3.51 


0.321 


1.17 


10.35 


88.5 


1.36 


3.33 


0.302 


1.11 


8.26 


74.4 


1.37 


3.15 


0.284 


1.04 


6.35 


61.1 


1.38 


2.96 


0.265 


0.97 


4.61 


47.5 


1.39 


2.78 


0.247 


0.90 


3.06 


34.0 


1.40 


2.60 


0.229 


0.84 


1.78 


21.2 



70 L , derived from UV data. In the case of a massive WD, the 
mass accretion rate is close to 1.1 ± 0.3 X 10~ 8 M yr^ 1 . The un- 
certainty in M is of the order of 23 percent and derives from the 
uncertainties in the distance, the correction for the inclination, 
and the value of Ri. 

An independent confirmation of these values for the mass ac- 
cretion rate can be obtained from the mass accretion rate versus 
boundary layer luminosity (M-Lbl) relation given in Table 2 of 
Patterson and Raymond (1985), if one assumes that Lbl ~ ^disk- 
Their Table 2 indicates that (for Mi = 1.0) a mass accretion rate 
M of 1.6 x 10~ 8 M yr~\ close to that derived above, would 
correspond to a luminosity of 1.3 x 10 35 erg s _1 . With the same 
M, a larger luminosity is expected by extrapolation to Mi val- 
ues close to 1.3 M , in close agreement with the value found for 
L disk (2.70 x 10 35 erg s" 1 ). 

The L1640-M relation reported in the same Table 2 does not, 
however, provide a consistent result. The average (de-reddened) 
flux on earth in the Hell 1640 line is 6.32 x 10~ 13 erg cm -2 
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s (Paper I) and the corresponding luminosity for d=3500 pc is 
9.3 x 10 32 erg s -1 . Table 2 of Patterson and Raymond (1985) 
indicates instead a Li 640 = 1.2 X 10 32 erg s _1 for M=10 18 g 
s~' = 1.587 x 10~ 8 M Q yr -1 , Mi=l, and extrapolation to higher 
Mi values does not appear to be able of reproducing the observed 
Hell luminosity. In T Pyx there is probably a nebular contribu- 
tion to L1640 that is added to that produced by the BL. It should 
also be emphasized that the Hell 1640 emission line exhibits 
significant changes, while the UV continuum remains constant 
in intensity and slope (see Paper I). This behavior casts doubts 
about the suitability of the L^q-M method for the specific case 
of T Pyx. 

Another estimate of M can be obtained by direct application 
of the M v -M relations reported in the literature. WLTO (1987) 
derived a relation between the absolute optical magnitude and 
the accretion rate for the "reference" disk model (Warner, 1987), 
after correction for the disk inclination: 

M corr. = _g M _ 5/3 . l og(Ml . M) ( 14 ) 

where Mi is in solar masses and M is in solar masses per year. In 
the case of T Pyx, for Mi in the range 1.25-1.4 M , one obtains 
values of M close to 5.0 x 10~ 8 M yr -1 . It should be noted, 
however, that in the M v -M relation given by WLTO there is 
no explicit dependence on Ri, which strongly depends on Mi, 
especially for massive WDs. If M is corrected for the factor 
0=(Ri/Mi)/(Ri ? /Mi o ) one obtains that M(1.25) = 3.50 x 10~ 8 
M yr" 1 and M(1.40) = 1.60 x 10~ 8 M yr" 1 . If one instead 
adopts the Paczynski and Schvarzenberg-Czerny (1980) correc- 
tion for i, the corresponding M values are 2.85 x 10~ 8 M Q yr~' 
and 1.22 x 10~ 8 M yr~' respectively. 

Lipkin et al. (2001) improved the M-M v relations presented 
by Retter and Leibowitz (1998) and Retter and Naylor (2000) 
and derived the new relation 

M„ = M- 4/3 x i ( 5 - 69 - - 4 ^) (15) 

where M 17 is the mass transfer rate in 10 17 g s _1 and M c l orr is the 
inclination-corrected absolute magnitude of the disk. The factor 
M[ 4/3 comes from the factor Ri/Mi in the assumption, generally 
valid for WDs with M < 1 M Q , that MR 3 =const. However, this 
Ri-Mj relation (polytropes) is not applicable if the WD is mas- 
sive because in this case Ri decreases with a far steeper slope 
as Mi increases (see Sect. 6 and Fig. 2). The adoption of the 
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standard Ri ~ Mj relation for massive WD has the effect of 
overestimating Ri and in turn overestimating M. This is espe- 
cially true for Mi ~ 1 .4, for which the difference between the Ri 
values estimated from the Ri ~ Mj relation and the more gen- 
eral Ri - Mi relations (Hamada and Salpeter 1961; Nauenberg 
1972; Anderson 1988; Politano et al. 1990; Livio 1994) reaches 
its maximum value (Fig. 2). 

The Lipkin et al. (2001) formula corrected for the proper 
Ri/Mi ratio in massive white dwarfs according to the <p fac- 
tors of Table 3, provides values of M = 1.38 x 10~ 8 and 0.59 
x 10~ 8 for white dwarfs with Mi= 1.25 M and Mi = 1.40 M , 
respectively. These values are about one half of the correspond- 
ing values obtained from the WLTO formula. It should be noted, 
however, that, in the derivation of Lipkin et al. (2001) formula 
there is the assumption that L v /Lfc ; ~ 0.14. In the case of T Pyx, 
this L v /Lfc ; ratio is lower by a factor of about two and therefore 
the M derived using their formula is underestimated by a similar 
amount. In conclusion, the various methods provide quite con- 
sistent results and we can confidently assume that the mass ac- 
cretion rate derived from "recent" UV and optical data is 1.86 ± 



Table 4. The basic quantities 



Quantity 


value 






-6.81 ±0.16 






6.7 ±0.1 




A v 


0.79 ± 0.06 




distance 


3500 ± 350 pc 




m v 


15.30 ± 0.05 




M v 
( 


1.79 ± 0.21 
25 ± 5 degrees 




Am v (i) 


0.74 ± 0.06 




M co rr 


2.53 ± 0.23 




Lw 


74 ± 15 L G 






70 ± 15 L Q 




M, 


~ 1.36 Mo 




M 


1.1 ±0.25 x 10" 


8 M yr 1 


M pre -OB 


2.2 ± 0.5 x 10- 8 


M yr 1 



0.4 x 10~ 8 M yr- 1 forM^l.25, and 0.84 + 0.2 x 10~ 8 M e y r - 1 
for Mi = 1.40 M . For an intermediate value of Mi = 1.33 M we 
obtain that M ~ 1.3 x 10~ 8 M yr~ l . 

Table 4 summarizes the adopted values of the basic parame- 
ters and their estimated errors. 

7. The pre-1966-outburst mass accretion rate 

Schaefer (2005) published a large collection of quiescent B- 
magnitudes for the recurrent novae T Pyx and U Sco using in- 
formation from archival plates, from the literature, and from his 
own collection of CCD magnitudes. The photometric data of T 
Pyx start in 1890 and therefore cover all of the known inter- 
outburst intervals. 

According to these data, during the inter-outburst phase in 
the years 1945-1966, T Pyx was a factor of 2 brighter in the B 
band than during the present quiescent phase after the 1967 out- 
burst. Since no significant changes in the (B-V) color index were 
found during both of these quiescence phases, one can safely de- 
duce that in 1945-1966 T Pyx was a factor of 2 brighter also in 
terms of visual and bolometric flux. Using the WLTO and the 
Lipkin et al. (2001) relations, one finds that the mass accretion 
rate in epochs pre-1966-outburst, hereinafter M pre -oB, was about 
twice the values of M for post 1967, obtained in the previous 
section. 

Therefore, we estimate that the mass accretion rate M pre -oB 
during the last pre-outburst interval (1945-1966) was between 
1.68 + 0.4 x 10~ 8 M yr 1 and 3.72 ± 0.8 x 10~ 8 M yr~ l , for 
M i = 1 .25 and M i = 1 .40 respectively. 

These are the values that are compared in Sects. 8 and 9 
with the theoretical models of nova. The post- 1967 accretion 
values (Table 3) are considered in Sect. 16 in the context of the 
"missing" outburst and the next expected outburst. 

8. The theoretical ignition mass and the accreted 
mass 

A nova outburst occurs when, due to the gradual accumulation 
of H-rich material on the surface of the white dwarf, the pres- 
sure at the bottom of the accreted layer becomes sufficiently high 
for nuclear ignition of H to begin ( Shara 1981, Fujimoto 1982, 
MacDonald 1983). Since the radius Ri of the white dwarf varies 
approximately as M[ I/3 for Mi < 1.0 M , the critical pressure 
for ignition 

P ign = (G ■ Mi ■ M ign )/(4nRt) (16) 
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corresponds to a critical ignition mass M igM that decreases ap- 

proximately as M, , while, for more massive WDs M, g „, de- 
creases with a steeper slope (see Table 3 and Fig. 2). In any case, 
massive white dwarfs need to accrete a small amount of mass to 
reach the critical conditions. 

For a given Mi value, the critical ignition mass has been cal- 
culated by various authors and the reported values (for a given 
Mi) differ from each other mostly in the choice of the critical 
pressure at the base of the accreted envelope, which varies be- 
tween 2.0 x 10 19 and 6.0 x 10 19 dynes crrT 2 (see Gehrz et al. 
1998; Stan-field et al. 1998; Livio and Truran 1992; Hernanz and 
Jose' 1998; Truran 1998). These studies indicate that a lower 
limit to Mj g „ (for a massive white dwarf with M\ close to 1 .4 M ) 
is in the range 2.0 - 4.0 x 10~ 6 M . However, as first pointed out 
by MacDonald (1983) and in more detail by Shara (1989) and 
Prialnik and Kovetz (1995), the behavior of a CN eruption and 
in particular the critical mass depends (apart from the WD mass) 
also on the mass accretion rate (since the WD is heated by ac- 
cretion) and the temperature of the isothermal white dwarf core. 

Townsley and Bildsten (2004) confirmed that the earlier pre- 
scriptions for ignition, based on the simple scaling M,- g „ oc R 4 
M^ 1 for a unique P, f „, are inadequate and that a system of a 
given Mi mass can have value of M, g „ that varies by a factor of 
10 for different M. At the high M values typical of most CVs, 
the critical pressure can decrease to values as low as 3 x 10 18 
dyn cm 2 , and the critical mass decreases accordingly. Thus, in 
a massive white dwarf accreting at significantly high rates, one 
expects a value for M Ig „ as low as 2-4 x 10~ 7 M G . 

The critical envelope mass Mt gn as a function of Mi and M 
can be numerically approximated to be : 

logM ign = -2.862 + 1.542 ■ M 1 " 1 436 Z«(1.429 - M\) + 

-0.l9(logM+ 10) L484 (17) 

(Kahabka and van Den Heuvel 2006) where Mi is in M and M 
is in M yr~ l . 

Table 5 indicates (for various Mi and hence Ri values) 
Mp re -oB, the theoretical Mj gn , the accreted mass M flccr =At- 
M l)re -oB, where At=22 yrs is the pre-1966 inter-outburst inter- 
val, and T=Mi g „/M pre -oB, that is, the expected recurrence time 
in years. The mass accretion rate and M, g „ were calculated using 
the quadratic fit for Ri as a function of Mi for massive white 
dwarfs derived in Sect. 6 (see also Table 3), and not the approxi- 

— 1/3 

mation Rj oc M, . Table 5 clearly shows that, after allowances 
for errors in the estimate of M pre -oB, the expected recurrence 
time r = Mig„J 'M pre -oB is close to the observed value (22 yr) for 
Mi ~ 1.36 - 1.38 M corresponding to M, g „ and M flccr in the 
range 3.0 to 6.0 x 10~ 7 M Q . 

This agrees with the estimate of Kato (1990) that the white 
dwarf mass for T Pyx is between 1.3 and 1.4 M Q , while in a 
subsequent paper Kato and Hachisu (1991) assumed Mi = 1.33 
M . 

We confirmed the results for M, f „ using the approximate re- 
lation for the ignition mass as a function of Mi and M given by 
Kolbetal. (2001): 

M ign = 4.4 x 10~ 4 • R\ ■ M- 1 ■ AT 1/3 (18) 

(where M,- g „ and Mi are in M , Rj is in 10 9 cm, and M is in 10~ 9 
M yr -1 ) and obtained M ig „ = 2.72 x 10~ 7 M . A similar result 
is also obtained by graphical interpolation in Fig. 5 of Kahabka 
and Van Den Heuvel (1997), which for M, ~ 1.37 and M ~ 2 x 



Table 5. Mi, the estimated pre-1967-outburst accretion rate 
Mpm-oB (for L</,' S/ t=140 L ), the theoretical ignition mass Mj gn , 
the accreted mass M flccr =22- M pre -oB and the expected recur- 
rence time T=MjgJ M pre -oB- 



Mi 


m pre-OB 




TVT 

"*ctccr 


T 


(M ) 


/ I a— 8 A If . T „— 1 \ 

(10 M Q yr ) 


{ 1 r\—7~\ /f \ 

(10 Mq) 


(10 Mq) 


(yrs) 


1.00 


7.32 


77.2 


16.10 


105.5 


1.05 


6.56 


67.6 


14.43 


103.0 


1.10 


5.84 


56.2 


12.84 


96.2 


1.15 


5.12 


45.7 


11.26 


89.2 


1.20 


4.40 


36.3 


9.68 


82.5 


1.25 


3.72 


25.7 


8.18 


69.1 


1.30 


3.02 


16.2 


6.64 


53.6 


1.33 


2.62 


11.0 


5.76 


42.0 


1.35 


2.34 


7.69 


5.15 


32.8 


1.36 


2.22 


6.14 


4.88 


27.7 


1.37 


2.08 


4.73 


4.58 


22.7 


1.38 


1.94 


3.44 


4.27 


17.7 


1.39 


1.80 


2.29 


3.96 


12.7 


1.40 


1.68 


1.33 


3.69 


7.9 



10~ 8 M yr~' provides an ignition mass close to 3.0 x 10~ 7 M 
(we note, however, that these parameters fall within a region in 
which only weak flashes are expected). In addition, we derived 
Pign as a function of M and the WD temperature from Fig. 3 of 
Yaron et al. (2005), and found that P ign ~ 4.8 x 10 18 for M = 
2 x 1O~ 8 M yr -1 (we note that in this regime there is a weak 
dependence on the WD temperature). The insertion of this value 
of Pig n in equation[8]gives M ig „ = 3.8 x 10~ 7 M . 

The consistency of all these results confirms that in the case 
of T Pyx the ignition mass was close to 4.5 x 10~ 7 M . 

In the recurrent nova T Pyx, the theoretical ignition mass and 
observed accreted mass are in excellent agreement for a massive 
WD, which provides new, independent support of the TNR the- 
ory. Theoretical expectations have received only limited confir- 
mation because studies of the system parameters (in U Sco, T Cr 
B) have provided contradictory results. 

9. Comparison with the nova models of Yaron et al. 
(2005) 

The realization that three basic and independent parameters, the 
white dwarf mass Mi, the temperature of its isothermal core T c , 
and the mass transfer rate M, control the behavior of a CN erup- 
tion, and the improvements in computer power and codes, has 
enabled an increase in sophistication in simulating a nova out- 
burst. Prialnik and Kovetz (1995) presented an extended grid 
of multicycle nova evolution models that have been extensively 
used by researchers. Each observed nova characteristic (e.g. 
peak luminosity, recurrence time, duration of the high luminos- 
ity phase, outburst amplitude, mass of ejecta, average outflow 
velocity, etc.) can be reproduced by a particular combinations 
of values of Mi ,T C , and M. Following this earlier study, Yaron et 
al. (2005) extended and refined the resolution in the grid of mod- 
els, including a considerable number of new parameter combi- 
nations. The full grid covers the entire range of observed char- 
acteristics, even those of peculiar objects. 

By matching the observed characteristics of a particular nova 
with its theoretical counterpart, it is therefore possible to derive 
information about the mass and temperature of the white dwarf 
and its average accretion rate. Therefore, the grids in Tables 2 
and 3 of Yaron et al. (2005) can be used to determine a set of 
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Table 6. A comparison between Yaron et al. (2005) grids and 
observations. 



Param. 


\/Ti\/|.-kl A 

ivioaei /\ 


ivioaei-D 


ivioaei-L, 


ivioaei-ij 


LJDS. 


Mi 


1.36 


1.33 


1.30 


1.40 


1.35±0.05 


M 


3.0 10" 8 


5.0 10" 8 


7.0 10~ 8 


1.0 io~ 8 


2.4±0.6 10" 8 


M ig „ 


6.6 10- 7 


1.03 10~ 6 


1.40 10- 6 


2.0 10~ 7 


-5.3 10~ 7 


M ej 


6.4 10~ 7 


0.93 10- 6 


1.30 10" 6 


2.0 10~ 7 


(see Sect. 10.2) 


T 


20 


20.6 


19.8 


20.2 


22 


Ampl. 


7.6 


6.9 


6.6 


8.4 


8.0 


Vexp.max 


-1500 


-720 


-572 


-1760 


-2000 


t 3 ,v, S 


12 


32 


45 


~6 


90 



fundamental parameters that provides the closest agreement with 
the observed data. 

In the case of T Pyx, a comparison of the observed parame- 
ters for the eruption phase (the recurrence time t , the outburst 
amplitude A, the decay rate t3, the expansion velocity v exp ) with 
the corresponding values in the grid clearly indicates that the 
closest agreement corresponds to values that are intermediate 
in between 1.25 and 1.40 for Mi, and -8.0 and -7.0 for logM; 
at these high M values, the effects of the core temperature T c 
are minor. Since a denser grid is not available, we calculated 
intermediate values by graphical and non-linear interpolation 
in the published values. A satisfactory agreement with the ob- 
served data comes from a model (model A) with Mi = 1.36 and 
logM=7.5 (M=3.2 x 10~ 8 ), which provides: t ~ 20, v exp ~ - 
1200, A - 7.6, t 3 - 13, and M ign ~ 5.0 x 10~ 7 . Following our 
interpolation of the data, Yaron (2007) kindly provided us with 
results obtained by execution of the code for other intermedi- 
ate models, that is, WD masses between 1.30 and 1.33, logM 
between -7 and -8, and different WD temperatures. The closest 
agreement with the observed values originates in a model with 
Mi = 1.33, Ti= 10 7 , logM=-7.3 (model B), although the model 
with Mi = 1.30, Ti= 10 7 , and logM=-7.15 (model C) is also ac- 
ceptable. In Table 6, we also include, as an illustration, a model 
directly obtained from the published Tables 2 and 3, for Mi = 1.4 
M Q and M = 1.0 x 10~ 8 (model D). 

The value for M, g „ (about 7 x 10~ 7 M Q ) that corresponds to 
the model of Yaron etal. (2005) with Mi = l. 33 and log M = 7.5, 
is close to that (~ 6.0 x 10~ 7 M ) estimated by the observed 
average M pre - B (about 3.0 x 10~ 8 M yr _1 ) and the observed 
inter-outburst recurrence time of 22 yrs. 

For T Pyx, the values of two of the independent, theoretical 
parameters (Mi and M) have been already constrained to within 
a limited range by theoretical considerations (Mi) or inferred by 
direct observations (M). Since the model grids require similar 
values for Mi and M, the results appear to be consistent. An ad- 
ditional fine tuning of the model data could certainly improve the 
fit but the overall agreement between models and observations 
can be considered satisfactory: an important parameter such as 
the recurrence time is reproduced quite well and for the other 
basic quantities the discrepancy is not large. 

A parameter that deserves special comment is t3. It is not ex- 
plicitly listed in the theoretical grid but can be associated to t m /, 
and was included by Yaron (2007) in specific calculations. The 
observed t3 value is about three times longer than the expected 
value. Models with such long t3 would correspond to higher M (g „ 
and therefore to far more massive ejecta. Models with a longer 
t3 (i.e. Model-C in Table 6) would be of too low v exp ,max an d too 
high M ign . 

In Sect. 10.2 we demonstrate that the spectroscopic observa- 
tions acquired during the last outburst in 1966 indicate the ejec- 



tion of a more massive shell than that expected from the simple 
estimate based on the quiescence M pre -OB and a times interval 
At of 22 years. 

Contini and Prialnik (1997) presented evolutionary calcula- 
tions of TNRs on the surface of an accreting WD that simu- 
lates the observed outburst characteristics of T Pyx and found 
that the model that reproduces the observed characteristics of 
T Pyx most accurately is obtained for an accretion rate of 10~ 7 
M Q yr~ l onto a 1.25 M WD; the (theoretical) mass ejected in a 
single outburst is expected to be close to 1.6 x 10~ 6 M . They 
also modeled the nebular spectrum and found that the theoreti- 
cal fluxes agree with those derived from observations (Williams 
1982) if a distance of - 2 kpc is assumed. However, our detailed 
study of the T Pyx luminosity and mass accretion rate suggested 
an M pre -oB value close to 3.7 x 10~ 8 M yr~\ onto a 1.25 M 
WD for d - 3.5 kpc. If d ~ 2.0 kpc, this would imply that M ~ 
1.2 x 10~ 8 , a factor of 10 lower than assumed by Contini and 
Prialnik (1997). We note that a distance as large as about 7.0 kpc 
is required to derive a luminosity correspomding to M ~ 10 -7 . 
See Sect. 1 1 . 1 for additional considerations of the distance. 



10. The photometric and spectroscopic history of 
the outbursts 

10.1. The long lasting optically thick phase 

The five recorded outbursts of T Pyx occurred in 1890, 1902, 
1920, 1944, and 1966, with a mean recurrence time of 19+5.3 
yrs (see WLTO, 1987). A common feature of these outbursts was 
the far longer optical decline time (t 3 = 90^) compared with that 
of other recurrent novae, which, with the exception of CI Aql (t3 
= 33 d) and IM Nor (t3 = 50 d), are generally much faster, with 
t3 of the order of days. 

The well studied outburst of Dec. 1966 exhibited a sharp ini- 
tial rise to a shoulder (pre-maximum halt on Dec. 10, 1966) in V 
of about 7.9, a near flat maximum close to V = 7.5 (which lasted 
about 30 days), a sharp peak at V=6.8 on Jan. 9, 1967, and a 
slow decline with t3 ~ 90 d . The first spectroscopic observations 
were obtained by Catchpole (1969), 12.6 days after the initial 
halt. The spectrum was characterized by sharp P Cyg features in 
the hydrogen lines and by other weak emission lines. The pres- 
ence of the absorption features endured until March 5, 1967, and 
since then they became quite "diffuse". Catchpole (1969), apart 
from the value of 850 kms~' at some epochs, reported outflow 
velocities of the order of -2000 kms -1 for both absorption and 
emission (permitted and forbidden) components. 

There is an observational gap of about 40 days in the spec- 
tra obtained by Catchpole, which was fortunately covered in part 
(from 09 Jan. 1967 to 16 Feb. 1967) by the spectroscopic data 
by Chincarini and Rosino (1969). They described the absorp- 
tion system in the Balmer lines (up to H12) as being particularly 
sharp and strong and noted an increase in the expansion velocity 
from -1535 km s _1 on Jan. 31, 1967, to -1760 km s _1 on Febr. 2, 
1967 and -1820 km s _1 on Febr. 6 1967. According to Chincarini 
and Rosino (1969), the mean radial velocity during the period 
Jan.31-Feb. 6 determined from all measurable absorption lines 
was -1810 + 40^8-'. 

Outflow velocities of the order of 1500 - 2200 kms -1 were 
also observed in previous outbursts. Adams and Joy (1920) re- 
ported the presence of "dark" components of radial velocity up 
to -2100 km s _1 a few days after maximum. Joy (1945) observed 
expansion velocities of 1700 kms 1 in forbidden emission lines 
in spectra taken some months after outburst. Similar values were 
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also found by Herbig (1945) and reported by Payne-Gaposchkin 
(1957). 

The "principal absorption system" is generally associated 
with the bulk of the mass ejected during outburst, and for most 
novae has a velocity close to that of the "nebular system" as 
determined from the width of the nebular emission lines (Payne- 
Gaposchkin 1957, Mc Laughlin 1956, Pottasch 1959). The v„ p 
value reported by Joy (1945) referred to observations during the 
nebular phase, three-four months after outburst. That the emis- 
sion lines in the nebular spectrum showed velocities of the same 
order as those deduced from the absorption lines of the principal 
spectrum is an indication of near constant expansion. 

We note that in the literature on T Pyx little attention has 
been paid to the fact that for almost three months after the ini- 
tial halt T Pyx showed a strong continuum with the presence 
of emission and absorption lines of similar strength. The persis- 
tence for almost three months (about t3) of displaced absorption 
components in the H (and Fell) lines, which are observable by 
combining the spectroscopic observations of Catchpole (1969) 
and Chincarini and Rosino (1969) for the 1966 outburst, indi- 
cates an optically thick phase of similar duration. 

In this respect, we recall that, before the outburst of 1966- 
1967, Mc Laughlin (1965) noted that T Pyx was exceptional 
among RNe since its photometric and spectroscopic behavior 
closely resembled that of a typical nova both close to maximum 
(where it remained for several weeks to within about a magni- 
tude) and in the nebular stage. 

10.2. The mass of the shell ejected in the optically thick 
phase 

The similarity between the spectroscopic and photometric char- 
acteristics of the outbursts of T Pyx and those of CNe, which 
allegedly eject about 10~ 4 - 10~ 5 M©, suggests in itself that dur- 
ing outburst T Pyx expelled a shell of comparable mass. 

Classical novae undergo an optically thick phase during 
which they resemble each other, a fact that can be explained by 
the same mechanism (i.e. flux redistribution) producing the spec- 
trophotometric light curve (Shore 1998, 2008). To achieve flux 
redistribution, the material must reach column densities of the 
order of 10 23 -10 24 cm , which corresponds to masses of about 
10~ 4 - 10~ 5 M . 

An optically thick stage also characterized the outbursts of T 
Pyx, as can be directly inferred from the lengthy period of time 
during which the optical magnitude was close to its maximum 
value, with t3 ~ 90 , and from the presence of absorption lines 
of HI and Fell, which lasted for at least 80 days. 

From the duration of the optically thick phase (associated 
to t2) and the observed v ex p, using simple assumptions, we can 
estimate the mass of the shell ejected during outburst. The outer 
radius of the shell can be estimated from the observed expansion 
velocity (y exp ~ 1 500 km s~ 1 ) and the time elapsed from outburst, 
assuming continuous ejection. This assumption is justified by 
the persistence of displaced absorption components with similar 
equivalent widths. The shell radius is: 

R ej ~ 7.7 x 10 14 cm ~ 1.1 X 1O 4 # . 

The corresponding shell volume is V~1.84 x 10 45 cm 3 . These 
values are conservative because the terminal velocity was prob- 
ably higher. To reproduce the optically thick stage recognized 
from the presence and persistence of the absorption lines, the 
column density must be of the order of 10 23 cm 2 . Therefore, 
the average density in the shell must be close to 10 8 - 10 9 cm" 3 , 



(we note that the density should scale as Rr 3 ). This density value 
agrees with the spectroscopic behavior and the presence of per- 
mitted emission lines only. If we assume that the ejecta are ho- 
mogeneous and consist of ionized hydrogen, the mass might be 
estimated as 

M H = N e m H V ~ 3.1 x I0 29 g = 1.5 x 1O~ 4 M . (19) 

This value is probably an upper limit because of the assumptions 
of continuous ejection and homogeneous shell. 

In an alternative approach, following Williams (1994), one 
can estimate the hydrogen column density produced by an ex- 
panding shell of mass 1.0 x 10~ 4 M : 

N H R = 3.0 X 10 52 • R- 2 [cm- 2 ]. (20) 

For T Pyx at day 60 (= t 2 ), we find that R 2 = 5.8 x 10 29 [cm 2 ], and 
N H R ~ 5.2 x 10 22 [crrT 2 ]. A mass of the ejecta higher than 10~ 4 
is therefore required to produce an optically thick stage until day 
60. 

It is well established that in novae, the mass of the ejected 
envelope is directly correlated with the optical decay time t 2 or t3 
(Livio, 1994). Therefore, an independent estimate of the ejected 
shell mass can be derived from the relation: 

logM ej = 0.274 + 0.197 • logt 2 - 4.355 + 0.283 (21) 

(Delia Valle et al. 2002). Even after considering the large un- 
certainties in this relationship, a t2=62 fl ' implies a mass for the 
envelope M e j ~ 10~ 4 M , which is similar to the mass ejected 
by classical novae. 

The data leading to Eq. [21] suffer from a large scatter. We 
suspect that the ejecta expansion velocity plays also an important 
role and should be included in the relation. Shore (2002, 2008) 
suggested an approximate scaling relation for the optically thick 
stage: 

M ej ~ 6.0 x 10~ 7 e N H ,24 V\ t\ M , (22) 

where V3 is the outflow velocity in 10 3 km s _1 , and e is the filling 
factor that, for this stage, can be assumed to be of the order of 
0.1. For t 3 =90, f=0.1, N H ,24 =0.1-1.0, and V 3 = 1.5, the derived 
values for M eJ are in the range 1.5 x 10~ 4 - 1.5 x 10~ 3 M Q . 

Finally, the ejected mass can be estimated using the follow- 
ing scaling law from Cassatella et al. (2005), which depends on 
the approximate assumption that the filling factor in novae ejecta 
is the same as in VI 668 Cyg: 

M ej ~ 0.044 ■ M Q /v exp , (23) 

where v exp is in km s _1 and the constant is set to the values 
of the ejected mass and the expansion velocity of VI 668 Cyg 
(Stickland et al. 1981). Using v ev; ,=1500 km s -1 , one finds that 
for T Pyx, M ej ~ 2.93 x 10~ 5 M . 

We recall that Kato and Hachisu ( 1 99 1 ), from their models of 
steady-state winds for a nova with Mi = 1.33 X=0.5 and Z=0.02 
and the observed t3, suggested the ejection of a massive envelope 
with M e j of about 10~ 5 M in a single outburst. 

All the previous results agree with the considerations of 
Shore (1998), who pointed out that, for a typical ejection ve- 
locity of about -2000 kms" 1 , a nova with an optical decline time 
longer than a week must eject a mass higher than 10~ 5 M . 

We also recall the fact that in classical novae close to maxi- 
mum, the Balmer lines develop P-Cygni profiles, which is a clear 
indication that a significant amount of material was ejected dur- 
ing the outburst (Starrfield, 1993). 

Therefore, all quantitative methods and the qualitative con- 
sideration of the photometric and spectral behavior of T Pyx dur- 
ing the outbursts indicate the presence of a massive envelope 
withM e ,- ~ 10~ 4 -10~ 5 M . 
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1 1 . The discrepancy between the mass of the thick 
shell and the ignition mass 

The ejection of a massive shell during the early (optically thick) 
outburst phase contrast significantly with the results of the UV 
+ optical observations during quiescence and the theoretical re- 
quirements for M ign , which imply a M pre - B ~ 2.2 x 1(T 8 and a 
total mass for the accreted shell M flccr of about ~ 5.0 x 1CT 7 M 
(see Sect. 8). It is also in contrast with the conclusions of Sect. 9 
that indicated that the closest agreement between the grid models 
and the observed properties of the system during outburst and Q 
corresponds to a model with M ~ 3.0 x 10~ 8 M yr _1 and M, g „ < 

I. 3 x 10~ 6 M (Table 6). During outburst, T Pyx has apparently 
ejected far more material than it has accreted. 

Studies of classical novae containing a massive WD indi- 
cated that these objects eject apparently more material than theo- 
retically predicted (Starrfield et al. 1998a; Starrfield et al. 1998b; 
Vanlandingham et al. 1996; Shore, 1998). These authors empha- 
sized the significant discrepancy between the observed mass of 
the ejecta and the predicted critical mass of accreted nova en- 
velopes for massive WDs (Mwd ^ 1-25 M ), the mass of the 
observed shell being one order of magnitude (or more) higher 
than that predicted by the models. 

For these CNe one could attribute the discrepancy to some 
inadequacy in the TNR models or in the methods to determine 
the nebular mass. In the case of T Pyx, the situation, however, 
differs because there is a serious mismatch between the shell 
mass indicated by the optical observations during outburst (M e j 
~ 10~ 4 10~ 5 M ) and that determined by the UV and optical 
observations during quiescence (which give M pre -OB ~ 2.2 x 
10~ 8 M yr~' and therefore M accr ~ M ign = 5.0 x 10~ 7 M ). 
Therefore, the mass ejected during outburst is about a factor of 
100 higher than both the theoretical ignition mass M ign and the 
total mass accreted before outburst, M accr . 

We note that using the post 1966 value of M instead of 
that inferred for the pre-outburst interval would produce an even 
larger discrepancy. 

II. 1. The role of the distance 

Although M e j does not depend on distance, M«, does (due to its 
dependence on both L disk and M pre _ 0B ), and therefore the mis- 
match between M e j and M accr depends crucially on the assumed 
distance. The uncertainty in the adopted distance was found in 
Sect. 2 to be of the order of 10 percent. However, we note, that 
even the adoption of an unlikely distance of, say, 10,000 pc (at 
about 20 cr from the estimated value) would only partially al- 
leviate this inconsistency, and at the expense of an uncomfort- 
ably high mass accretion rate, well within the range for the onset 
of steady burning. This would imply characteristic temperatures 
and luminosities that are not observed (see also Sect. 14 for a fur- 
ther discussion). A larger distance would also necessarily imply 
that T Pyx was super-Eddington at maximum, a circumstance 
that appears unlikely due to its slow photometric and spectro- 
scopic developments during outburst. 

A lower distance would correspond to a lower values of Lj, ^, 
and hence M pre -oB and M flccr , to values that are theoretically in- 
compatible with the occurrence of outbursts with an average in- 
terval of 22 years. It would also exacerbate the discrepancy be- 
tween the low value of M accr obtained from the UV observations 
(and the models) and the apparently high mass of the ejecta, as 
inferred from the behavior during outburst. Therefore, there is 



not much leeway to invoke a different distance to explain the 
discordance. 



12. The nebula revisited 

The nebula surrounding T Pyx has been the target of several 
spectroscopic and imaging observations. Duerbeck und Seitter 
(1979) first reported the presence of a strong nebulosity around 
T Pyx, with radius r ~ 5", whose origin was tentatively attributed 
to the 1966 outburst and whose strength was described as un- 
usual. By the assumption of an outburst expansion velocity of 
-900 kms -1 , a low distance (~ 600 pes) was derived. 

Williams (1982) obtained spectral scans of the northern por- 
tion of the nebula of T Pyx. The spectrum was similar to that of 
as a typical PN, probably photoionized by radiation from the hot 
remnant, and lacked the strong CNO enhancements characteris- 
tic of the ejecta of classical novae. 

Comparing images acquired in 1979.0 and 1982.9, Seitter (1987) 
found that the nebula did not increase in size during that time 
interval. 

Shara et al. (1989) from deep narrowband CCD images con- 
firmed the faint, extended H ff + Nil halo (twice as large as the in- 
ner nebula), first reported by Duerbeck (1987). A smooth, small 
[OIII] nebula with r ~ 2" was also found. Shara et al. (1989) also 
compared the relative sizes of the main nebula with r ~ 5" in 
1985 and in 1979 but failed as well to find any detectable expan- 
sion during the 6 year interval, confirming the finding of Seitter 
(1987). High resolution imagery data from HST was obtained 
by Shara et al. (1997). The nebula was resolved into more than 
two thousand individual knots, and a comparison between im- 
ages taken at four epochs indicated that these individual knots 
retained a similar pattern, without any evidence of expansion. 
These data confirm the apparent stationarity in the 10" diameter 
nebula suggested by previous observations. Shara et al. (1997) 
found an upper limit of 40-(1500/d) (km s _1 ) for the expansion 
velocity of the knots. They also detected nine distinct peaks in 
the brightness distribution, an indication that a multiple nebula 
model was required. 

Many studies, disappointingly inconclusive, addressed the 
problem of the mass of the nebula. In this respect, it should 
be noted that large uncertainties are generally associated with 
estimates of the mass of novae ejecta, since the mass estimate 
depends critically on quantities that are not reliably measured, 
e.g.: distance, electron density, ionization structure in the nebula, 
geometry, and filling factor. Therefore, it is unsurprising that a 
range of masses for the nebula of T Pyx has been proposed in 
the literature. 

From the nebular H^ intensity measured by Williams (1982), 
WLTO (1987) obtained a lower limit to M. neb of 10~ 6 M , while 
Shara (1989), using the H^ intensity with the requirement e <1 
for the filling factor derived an upper limit of 1.0 x 10~ 4 M . 
From the intensity of the U a and [Nil] lines, Seitter (1987) found 
a mass close to 8.0 x 10~ 5 M . From the H a flux and consider- 
ations based on HST imagery, Shara et al. (1997) obtained 1.3 
X 10~ 6 M to be the most reliable estimate for the nebula mass 
(with an assumed distance of 1500 pc), the electron density be- 
ing, allegedly, the main uncertainty factor. 

We add one more estimate for the mass of the nebula, based 
on the H B flux obtained by Williams (1982) scaled to the entire 
nebula, after correction for the new reddening and distance. We 
obtain : 

F B ~ 4.24 x 10~ 14 erg cirT 2 s" 1 and L B = 6.22 x 10 31 erg s _1 . 
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with a relative error of about 25 percent, arising from the un- 
certainties in the distance and the flux. By combining the two 
common relations: 

L p = 1.24 X 1(T 25 N e ■ N + ■ (4/3) n R 3 neb ■ e (24) 
and 

M neh = M -N + - m H ■ (4/3) 7T R 3 neb ■ e, (25) 

(where //=1.4 is the mean atomic weight and e is the filling fac- 
tor), one finds that the mass of the nebula (independent of e and 

Rneb) IS 

U 0.584 

M neb = 18.67 ^ g = M . (26) 

N c is inaccurately determined because, as noted by Williams 
(1982), the spectrum is too poor in number of lines to constrain 
the physical parameters of the nebula accurately. However, from 
the presence of a quite strong emission feature near 3722 A, at- 
tributable to [Oil] 3726.1 and 3728.8 one can estimate that the 
nebula density is less than 3.0 x 10 3 cirT 3 , since these two for- 
bidden lines have critical density log N c e r " = 3.5 and 2.8, respec- 
tively. If we assume that N c ~ 10 3 otT 3 , a value also adopted by 
Shara et al. (1997) from other considerations, we find that M ne i, 
~ 5.84 x 10~ 4 M Q . 

A value of this order for the total mass of the nebula, which 
has apparently increased from the contribution of several suc- 
cessive outbursts, agrees with the ejection during outburst of a 
massive shell (~ 10 -5 M ) as suggested in Sect. 10.2. 

The fact that the nebula of T Pyx is still clearly observed 
despite its large distance seems hardly compatible with a mass 
of 10~ 6 M derived in previous studies. In the most well studied 
CNe, which on average are closer in distance than T Pyx, the 
ejected nebula has, at best, a similar strength and/or is barely 
evident a few years after outburst. We think that the peculiar 
strength of the nebula of T Pyx may be explained by the fact that 
most of the gas ejected in successive outbursts has accumulated 
into a nearly stationary envelope that is strongly irradiated by a 
more than average luminous UV central source. 

Alternatively, one can guess that the observed nebula of T 
Pyx was produced in a peculiar event and is not associated with 
the recorded and/or previous outbursts. Its apparent stationarity 
and the lack of changes over a timescale of about fifteen years 
(Shara et al., 1997) supports this interpretation. Williams (1982) 
noted that the spectrum of the nebula around T Pyx is almost 
similar to typical planetary nebula, of approximately solar com- 
position. 

We emphasize that the optically thick shell discussed in 
Sects. 10.1 and 10.2 was observed only spectroscopically dur- 
ing the outburst phases, and that there is no definite, direct link 
with the extended nebula. Simple calculations indicate that, if 
d=3500 pc, the angular radius of the shell, after ten years of 
constant expansion at 1500 km s" 1 , would be less than 1". 

12.1. The filling factor 

If the radius of the nebula R ne b is known, from the observed an- 
gular radius r (< 10") and the distance, one can obtain the filling 
factor e with the help of eq. [24] After insertion of the values for 
Lp = 6.22 x 10 31 erg s -1 , R„ eh ~ 5.2 x 10 17 cm, and N e ~ 1.1 
N + , one obtains that 
e ■ N 2 e ~ 1062.7. 



Therefore, for N e ~10 3 cirT 3 , e is close to 10~ 3 . A filling factor 
of ~ 10~ 3 can be obtained independently from the relation of 
Harrison and Stringfellow (1994): 

e ~ NdVJVneb = 8N d ll2 y V2 , (27) 

where y is the fraction of the spherical nebula intercepted by the 
clumps. From the figures of Shara et al. (1997), we estimate that 
y ~ 0.05 and, if the number of clumps N c ; is ~ 2000 (Shara et 
al. 1997), we find that e ~ 2 x 10 3 . We note that the filling 
factors estimated for the ejecta of other novae have a wide range 
of values, from 4 x 10~' (Vanlandingham et al., 1999 for Nova 
LMC 1990 no.l), to 10" 6 (Saizar and Ferland (1994) for nova 
QU Vul. Williams et al. (1981) found an intermediate value, e 
~ 2 x 10~ 3 for the recurrent nova U Sco. Mason et al. (2005) 
found e ~ lO^-lO 1 for nova SMC 2001 and nova LMC 2002, 
and Balman and Oegelman (1999) estimated e ~ 5 x 10~ 3 -3 x 
10 _1 for the shell of GK Per. Mason et al. (2005) have suggested 
e~ 10- 5 -10- 2 forTPyx. 

13. The mass balance and the SNIa connection 

It is accepted that Type la supernovae represent the complete 
thermonuclear disruption of mass-accreting white dwarfs that 
reach the Chandrasekhar limit by accretion (Nomoto et al. 1984; 
Woosley and Weaver 1986). Within this general framework, 
there exist single degenerate models (Whelan and Iben, 1973) 
in which a WD accretes from a non-degenerate companion, and 
double degenerate models (e.g. Iben and Tutukov, 1984) that in- 
volve the merger between two WDs. Hachisu (2002, 2003) and 
Hachisu and Kato (2001, 2002) proposed a unified picture of bi- 
nary evolution to SNe la in which recurrent novae could be un- 
derstood to be part of the evolutionary stages of supersoft X-ray 
sources and symbiotic channels to SNe la. 

The mass Mi must be close to the Chandrasekhar limit be- 
fore a recurrent nova can become a SN la, and the WD must 
also increase in the long term, after many cycles of accretion 
and ejection (Starrfield et al., 1985). However, in the case of T 
Pyx, even if Mi were close to the limit, which is not clearly es- 
tablished since Mi appears to be close to 1.37 M , the results of 
the previous sections indicate that the mass balance situation is 
unclear: 

On the one hand, the photometric and spectroscopic behav- 
ior close to outburst, as mentioned in Sect. 10.2, appear to be 
consistent with the ejection of a rather massive shell, while the 
UV data and theoretical models (at the specific Mi and M val- 
ues) indicate that the ignition mass is low (~5.0 x 10~ 7 M Q ). 
This indicates that during outburst the WD ejects more material 
than it accumulates and that a secular decrease in the mass of the 
white dwarf is expected. Therefore, evolution to become a SNIa 
appears to be excluded. 

On the other hand, the ejection of a more-massive-than- 
accreted shell is apparently in contrast with the observational 
evidence that the chemical composition of the T Pyx nebula is 
close to solar (Williams, 1982). This appears to exclude any ero- 
sion of the white dwarf and implies that the white dwarf does 
not lose mass after cycles of accretion and ejection. However, 
this presupposes that the chemical composition of the observed 
nebula is representative of that of a single shell ejected during 
outburst. 

It is unclear whether these substantial discrepancies origi- 
nates in flaws in the theoretical assumptions or the interpretation 
of the observations. They certainly highlight the need for accu- 
rate values of the most critical parameters of this recurrent nova, 
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i.e. the mass and chemical composition of the shell ejected dur- 
ing outburst. 

In any case, the behavior of T Pyx raises several doubts 
about the common assumption that in recurrent novae far less 
material is ejected during the outburst than is accreted by the 
white dwarf, and therefore that its mass increases toward the 
Chandrasekhar limit (Starrfield, 2002). In assessing whether re- 
current novae could be progenitors of Type la supernovae, we 
note that Hachisu and Kato (2001), discovered that, among six 
recurrent novae, only T Pyx is offset significantly from the re- 
gion occupied in an a orbital period - donor mass plane (their 
Fig. 3) by the progenitors of SNe la. 



14. Neither a supersoft X-ray source, nor assisted 
suicide 

To explain the alleged extremely blue color of T Pyx in quies- 
cence, WLTO (1987) proposed that nuclear burning continues 
even during its Q state, consistent with the slow outburst de- 
velopment, which suggests that the accreted envelope was only 
weakly degenerate at the onset of TNR. Patterson et al. (1998) at- 
tributed the luminosity of T Pyx (and V Sagittae) to quasi-steady 
thermonuclear burning and suggested the object to be included 
in the class of the supersoft X-ray sources. 

However, the color of T Pyx (B-V) = -0.26 used in 
these studies is based on a significant overcorrection for the 
reddening, assumed to be E B _y=0.36, instead of the correct 
value E fi _y=0.25 (see Paper I and earlier communications, e.g. 
Gilmozzi et al., 1998). It is unfortunate that both the value (B- 
V)„=-0.26 of Patterson et al. (1998) and the statement about the 
"extremely blue color" of T Pyx was adopted widely in the lit- 
erature (see for example Anupama 2002, and Parthasaraty et al. 
2007). We note, incidentally, that in the same paper Patterson et 
al. (1998) assumed too high a reddening correction (Efl_y=0.33) 
for V Sge; the IUE data suggest, instead a value close to 0.23. 
We recall that for T Pyx the observed (B-V) is about 0.14 + 
0.04 (WLTO 1987; Bruch and Engel 1994; Downes et al. 1997; 
Schaefer 2005; see also Table 7). This would imply that (B-V)„ ~ 
-0.1 1, close to the value (B-V) = -0.06 given by Szkody (1994). 

Patterson et al.(1998) assumed that M v =1.3 and after a sig- 
nificant bolometric correction (based on the assumption of an 
extremely hot source being present, a consequence of the over- 
estimate of the reddening), derived a quiescent bolometric lumi- 
nosity higher than 10 36 erg s _1 , which considered to be a true 
lower limit. This encouraged them to invoke nuclear burning on 
the surface of the WD as the main power source, considering the 
disturbingly high M ( > 10~ 7 M Q yr~') required in the case of 
pure accretion power. 

However, the IUE and the optical data do not appear to be re- 
produced by the model depicted by Patterson et al. (1998), since 
the following observational evidence contradicts with their con- 
clusions: 

1 . One of the main results of Paper I was that the de-reddened 
UVOIR continuum of T Pyx is reproduced well by a single 
power-law F A oc A a with a slope a = -2.33, representative 
of a steady accretion disk. Alternatively, as shown in Paper 
1, the UV continuum can be well fitted by a blackbody of 
temperature 34000 K. There is no way to reconcile these firm 
observational results with the presence of a supersoft source 
of typical T ~ 3.0 X 10 5 K, since its expected slope in the 
UV region (a = -3.80) would be inconsistent with the UV 
observations. 
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Fig. 3. Comparison between the reddening-corrected S WP spec- 
tra of the supersoft source V Sge (continuous line) and of T Pyx 
(dashed line). The continuum of T Pyx has been scaled to that 
of V Sge (multiplied by a factor 40). Note the prominence of the 
high ionization lines of NV A 1240, CIV A 1550, Hell A 1640, 
and NIV A 1719 in V Sge, in contrast with their moderate inten- 
sity or absence in T Pyx. 

2. As shown in Sect. 5, the IUE observations in 1980-1996 and 
the optical and IR photometric data indicate that Lj,-. s/ t=2.7 
x 10 35 erg s _1 , which corresponds to emission mainly in the 
directly observed UV range. Therefore, in the absence of any 
direct or indirect evidence of a hot source, it is unlikely that 
L disk is higher than 10 36 erg s" 1 . 

3. Under steady-state hydrogen-burning conditions, the accre- 
tion rate M steady can be estimated to be: 

M steady ~ 3.7 x 10- 7 (M! - O.4)M yr~\ (28) 

(Hachisu and Kato (2001), which is valid for hydrogen con- 
tent X=0.7 In the case of T Pyx, one calculates M steady ~ 3.2 
x 10~ 7 M G yr~ l . This value is almost a factor of 30 higher 
than the value of M ~ 1.1 x 10 -8 M yr~ l obtained from the 
IUE data for the post outburst phase. 

4. After correction for inclination, the apparent absolute mag- 
nitude of M v =1.79 corresponds to a 4-n averaged absolute 
magnitude of M(; orr =2.53, to be compared with the average 
value M™" ~ 4.0 for ex-novae (Warner, 1995). This implies 
that T Pyx is more luminous than other ex-novae, as men- 
tioned in Sect. 4, but is not "extremely bright". We can ex- 
plain why T Pyx is brighter than other novae in terms of 
repeated nova eruptions and heating of the primary, which 
triggers irradiation of the secondary and produces a higher 
than average M. 

5. The emission line spectrum of T Pyx is not that of an high 
excitation object: NV A 1240 is nearly absent, Hell A 1640 is 
weaker than CIV A 1550 and barely present in several spec- 
tra, the OIV lines close to A 1405 are absent. A comparison 
between the IUE spectra of T Pyx and V Sge (Fig. 3) clearly 
shows remarkable differences and a much lower excitation 
character in T Pyx. From an inspection of IUE spectra of 
several CVs, we also found that the spectrum of T Pyx is 
similar to that of the old novae V533 Her and V603 Aql and 
to some spectra of the intermediate polar TV Col. The old 
nova RR Pic definitely shows higher excitation than T Pyx, 
with far stronger NV and Hell emissions. 
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Other findings considerations also exclude the H-burning- 
bloated-WD hypothesis : 

1 . The complexity of the optical photometric behavior in T Pyx 
(Shaefer et al.,1992; Patterson et al., 1998) is difficult to ex- 
plain if the bulk of the luminosity originates in a spherically 
symmetric radiation source associated with H-burning on top 
of a bloated white dwarf. 

2. If the majority of the gas undergoes steady burning, it would 
be difficult to understand how the remainder that accumu- 
lates in the degenerate envelope could burn explosively every 
20 years. 

3. The outburst amplitude of T Pyx, ~ 8.0 magnitudes, is close 
to that found for classical novae of similar t3 and similar 
system inclination (Warner, 1995, 2008). If, as assumed by 
Patterson et al. (1998), the luminosity during quiescence is 
greater than 10 36 erg s" 1 then, with an outburst amplitude 
of about 8 magnitudes, the luminosity would reach 10 40 erg 
s _1 , implying that T Pyx is an object intermediate between 
a nova and a SN la. 

To support the hypothesis of steady nuclear burning 
Patterson et al. (1998) considered all CVs of comparable P or t, 
and deduced that the M in T Pyx was a factor of 5000 higher 
than in other CVs of similar P on j,. One should compare T Pyx 
with objects that are similar, that is, with recent novae. After 
outburst, the nova system remains in an excited state and M in- 
creases due to the irradiation of the secondary. Patterson et al. 
(1998) correctly excluded from their Fig. 14 all old novae within 
30 yr of the outburst because of their systematically too high 
luminosity levels. Adopting the same line of reasoning, T Pyx 
should also have been excluded. In this respect, we note that the 
M of T Pyx is only slightly higher than that observed in recent 
ex-novae (e.g. RR Pic, V8410ph, HR Del, etc) (Selvelli, 2004). 

Based on the conclusions of Patterson et al. (1998) of an the 
extremely high luminosity (L/, fl ; far higher than 10 36 erg s~'), 
Knigge et al (2000) investigated in detail the evolution of the T 
Pyx system and proposed that the system is a wind-driven super- 
soft X-ray source. In this scenario, a strong, radiation-induced 
wind is excited from the secondary star, and increases the rate of 
binary evolution, causing the system to destroy itself either by 
evaporation of the secondary star or in a Type la SN if the WD 
reaches the Chandrasekhar limit. Knigge et al. (2000) therefore 
proposed that either the primary, the secondary, or both stars may 
be committing assisted stellar suicide 

This scenario is, admittedly, highly speculative, and depends 
crucially on the unsubstantiated assumption that both the tem- 
perature and luminosity of T Pyx are extremely high. The IUE 
and optical data are instead consistent with a more conventional 
scenario of accretion power, as in other CVs, and we confidently 
predict that, fortunately, any form of suicide in the near future is 
extremely unlikely. 

Finally, we note that Greiner et al. (2003) did not find T Pyx 
to be a supersoft X-ray source, and that T Pyx does not appear 
in NASA's HEASARC tool (a master compilation of EUV and 
X-ray databases). 

15. The XMM observations 

While this work was close to completion, the data for X-ray ob- 
servations of T Pyx by XMM became publicly available. This 
prompted us to perform a preliminary analysis of the data to ver- 
ify the presence or absence of a supersoft source. 

T Pyx was observed by XMM-Newton on November 10 
2006. All the three EPIC cameras were operated in Full Frame 
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Fig. 4. The XMM-Newton EPIC-pn spectrum of T Pyx (bot- 
tom) compared with the simulations of a 20 ksec exposure of a 
blackbody of 2.4 x 10 5 K and a luminosity of 1 x 10 37 erg s" 1 
computed with two assumptions: a distance of 3500 pc and a 
reddening E B _y=0.25 (values assumed in this paper, dots), and a 
distance of 3000 pc and a reddening of 0.4 (values assumed by 
Knigge et al., 2000, continuous line). The three spectra shown 
here have been re-binned to 20 counts per bin. 



mode with the Medium filter. The total useful exposure time af- 
ter filtering for high radiation periods was 22.1 ksec. Optical 
Monitor data were taken simultaneously with the X-ray obser- 
vations. The values (see Table 7) are consistent with the values 
given in Paper I (IUE and optical observations) and confirm the 
stability of the SED with time. The reduction of the XMM EPIC 
data was carried out with SAS version 7.1, using standard meth- 
ods. T Pyx was detected as a faint source that had an observed 
EPIC-pn count-rate of 8.5 x 10~ 3 cts s _1 and emission over the 
complete range 0.2-8 keV. 

Figure 4 shows the XMM-Newton EPIC-pn spectrum of T Pyx 
compared to a blackbody of 2.4 x 10 5 K and a luminosity of 1.0 
x 10 37 erg s _1 . The blackbody spectrum was simulated assuming 
an exposure time of 20 ksec, similar to that of the data, with two 
assumptions: a distance of 3500 pc and a reddening E B _y=0.25 
(as in this paper), and a distance of 3000 pc and a reddening 
of 0.4 (as assumed by Knigge et al., 2000). The presence of a 
supersoft component with a temperature of 2.4 x 10 s K and a 
luminosity of 1.0 x 10 37 erg s _1 , whose existence was postulated 
by both Patterson et al. (1998) and Knigge et al.(2000), can be 
definitely excluded. Such a bright component would be easily 
visible at soft energies (i.e. below 0.5 keV), and this is not the 
case. Any supersoft emission, if present, would be several orders 
of magnitude fainter than expected. 

The data are compatible with the presence of a relatively hard 
source of low luminosity, but their low statistical quality (227 
counts in the range 0.2-8. keV) does not allows us to perform a 
reliable fit. We also note that extrapolation of the observed IUE 
spectrum (power-law) into the X-ray range would result in an 
unrealistically high flux, several orders of magnitude higher than 
observed. The existence of such a high flux was already excluded 
in Paper I on the basis of the intensity of the Hell 1640 A line. 
We therefore ask from which process or in which source the ob- 
served X-ray emission originates ? In a CV (nova or quiescent 
nova) X-rays may be produced by at least three different mecha- 
nisms or regions (see Kuulkers et al. 2006, and Krautter, 2008): 
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1. Thermal emission from the hot white dwarf nova remnant in 
the late outburst phases. The nova becomes a strong X-Ray 
emitter with a soft SED. 

2. Emission from the inner accretion disk (BL) or the accre- 
tion columns (in magnetic CVs). One expects to observe a 
typical X-ray emission of CVs in quiescence, with a thermal 
bremsstrahlung spectrum. 

3. Shocks in the circumstellar medium surrounding the nova 
system where the expanding nova shell and/or a nova wind 
interact with each other or with pre-existing CS material. 
The expected hard X-ray spectrum originates from ther- 
mal bremsstrahlung with kTeff temperatures 0.2-15 keV, and 
with Lx ~ 10 33 -10 34 erg s" 1 (O'Brien et al. 1994). 

4. The corona of an active dwarf M star companion. 

In T Pyx, as already mentioned, the XMM data would ex- 
clude case 1. corresponding to a strong SSS source. Case 2 can 
also be ruled out because, for a CV accreting at the quite high 
rates of T Pyx, the BL is optically thick and one would expect 
to observe (see Patterson and Raymond 1985) only soft X-ray 
emission with luminosity comparable to the disk luminosity, Lbi 
~ Ldisk ~ 2.7 x 10 35 erg s _1 . For case 4, we note that a high 
level of coronal heating at values of L x fLb i ~ 10~ 3 persists for 
the active dMe stars, and the emission was detected in surveys 
of nearby stars (see Giampapa and Fleming, 2002), although any 
detection would be impossible at the distance of T Pyx. 

We suggest that the most likely origin of the observed hard 
X-ray emission is from shocks within the circumstellar enve- 
lope. We note that in GK Per, Balman et al. (2006) detected 
hard X-ray emission by direct imaging with Chandra. The to- 
tal X-ray spectrum of the nebula consists of two thermal promi- 
nent components of emission. GK Per has a large amount of CS 
material, which is most likely a residual of a planetary nebula 
phase, and the shell remnant shows a clumpy structure similar 
to that observed in T Pyx by Shara et al. (1997) with HST. We 
recall that the studies by Orio et al. (2001) by Orio (2004) of 
the X-ray emission from classical and recurrent novae demon- 
strated that emission from shocked ejecta is expected to last 
about two years, but may last for up to a century, if, for example, 
there is pre-existing circumstellar material (as in the case of GK 
Per). Hernanz and Sala (2007) reported on X-ray observations of 
V4633 Sgr performed with XMM-Newton between 2.6 and 3.5 
yr after outburst. The X-ray spectrum is dominated by thermal 
plasma emission, which most probably originated in the shock 
heated ejecta. 

Unfortunately, the limited spatial resolution of the available 
XMM data do not enable any spatially-resolved study to be com- 
pleted, because the pixel size is about 4 arcsec compared with the 
optical radius of the nebula which is about 5 arcsec. 

Also the XMM observations, excluding the possibilities of 
continuous burning and the supersoft source scenario (a massive 
white dwarf accreting at high rates), appear to exclude T Pyx be- 
coming a SN la by means of the supersoft X-ray source channel 
described by Hachisu (2002, 2003). 



16. The recurrence time and the next, long-awaited 
outburst. 

As reported in Paper I, we started an observing program in 1986 
with IUE to monitor T Pyx prior to (and during) the expected 
next outburst that was supposed to occur in the late eighties of 
the last millennium. Unfortunately, the star successfully man- 
aged to postpone the long-awaited outburst, and at the present 



Table 7. XMM observations compared with IUE fluxes and 
ground based photometry 



Filter 


Xeff 

(A) 


Observed flux 


IUE flux 


magn. 




(1 x 10~ 14 ) 


(1 x 10~ 14 ) 




UVM2 


2310 


0.92 


0.80 




UVW1 


2910 


0.86 


0.92 




u 


3440 


0.70 




14.37 


B 


4500 


0.38 




15.59 


V 


5430 


0.24 




15.49 



time (2008) has surpassed by eighteen years the longest inter- 
outburst interval so far recorded (24 yrs). 

As mentioned in Sect. 7, Schaefer (2005) published the re- 
sults of a study of the inter-outburst interval in the recurrent no- 
vae T Pyx and U Sco. From an analysis of the available data, he 
found that the two novae are relatively bright during short inter- 
eruption intervals and dim during long intervals, suggesting that 
the product of the inter-eruption interval times the average bolo- 
metrically corrected flux is a constant. Therefore, in the case of 
T Pyx, the lack of the post- 1967 outburst is explained by a lower 
luminosity and therefore a lower M. From the decline in the ob- 
served quiescent B magnitude in the time intervals before and 
after the 1966 outburst, Schaefer (2005) also predicted that the 
next outburst of T Pyx will occur around 2052. 

With the help of considerations in Sect. 8 and the data in 
Tables 3 and 5, we can further investigate this prediction. The 
recurrence time can be estimated from the (theoretical) M ign and 
the observed mass accretion rate. M (g „ depends mainly on Mi 
and M, while M, for a given L^, is a function also of Mi and 
Ri, which, in turn, is a function of Mi. Tables 3 and 5 list M (g „, 
M, and the recurrence time t=M,- i? „/M (years) for various Mi 
values. Table 5 (pre-1967 M values) clearly shows that the ob- 
served inter-outburst interval (22 years) corresponds to Mi val- 
ues ~ 1.36-1.38 M Q . Table 3, which contains post 1967-outburst 
M values indicates that the observed interval, which, so far, is 
longer than 42 years, corresponds to Mi < 1.38 M Q . The most 
likely value of Mi is therefore close to 1.37 M Q . 

We note that a reduction by a factor two in the mass accretion 
rate corresponds to an increase in the expected recurrence time 
t by a larger factor because, for a given Mi, M, g „ increases as 
M decreases. For the relevant values of T Pyx, the decrease by 
a factor of 2 in M is accompanied by an increase by a factor of 
about 50 percent in M ign . For the next outburst, therefore one 
expects an increase in the inter-outburst interval t by a factor 
of approximately 3.0, to values near 60 years (see Table 3). Our 
prediction for the next outburst date is therefore around A.D. 
2025. With this new date, contrasted with that of Schaefer (A.D. 
2052), we (or at least some of us) feel a bit more confident about 
the chance of personally testing this prediction. 

However, given the uncertainties in M and M ign , the possibil- 
ity of a more imminent outburst cannot be ruled out. In this case, 
X-ray and other observations during the first outburst stages will 
be of paramount importance in determining the mass ejected in 
a single event. 

17. Summary and conclusions 

We have accurately determined, from UV and other observa- 
tions, the accretion disk luminosity of T Pyx during both the 
pre- and post-1966 inter-outburst phases. For Mi ~ 1.37 M Q , 
we have found that M pre -OB ~ 2.2 x 10~ 8 M yr -1 . By combin- 
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ing the measured accretion rate with the duration of the inter- 
outburst phase (At = 22 yrs), the total accreted mass is inferred 
to be M accr =M pre - B ■ At ~ 5.2 x 10~ 7 M . This value is in ex- 
cellent agreement with the theoretical ignition mass (M, g „) ~ 5.0 
x 1(T 7 M G expected for a massive white dwarf accreting at the 
quoted rate. Therefore, both the time interval between the last 
two outburst and the absence of the awaited post- 1967 outburst 
(due to the lower M in the post- 1967 time interval) are explained 
in a self-consistent way. 

This is the first reliable determination of the mass accreted 
prior to a nova outburst, M accr , owing to the dominance of the 
accretion disk luminosity over that of the secondary star at UV, 
optical and IR wavelengths, as well as good observational cov- 
erage during the inter-outburst phases. Unfortunately, M accr can- 
not be confidently determined in other cases, such as classical 
novae, because of their long inter-outburst interval, nor in other 
recurrent novae, due to the faintness of the source, the lack of 
systematic UV observations, or the dominance of light from the 
giant companion over that from the accretion disk. 

In T Pyx, the consistency between the observed M accr and 
the theoretical M, g „ supports the good quality of the observations 
and the reliability of the models and represents a new, direct con- 
firmation of the validity of the TNR theory, which associates a 
massive white dwarf with the recurrent nova phenomenon. A 
detailed comparison of the observed parameters with the the- 
oretical grids of Yaron et al. (2005) indicates that the closest 
agreement is obtained with a models of a rather massive white 
dwarf (Mi ~ 1.25-1.40 M Q ) that accretes at high M rates (M ~ 
10 8 -10 7 M G yr -1 ). However, no combination of the theoretical 
parameters can reproduce the observed values reliably, t3 being 
the most difficult parameter to describe. 

The literature data of the spectroscopic and photometric evo- 
lution during the outbursts of T Pyx clearly indicate the occur- 
rence of an optically thick phase that lasted about three months. 
This implies an ejected mass of M e j ~ 10~ 5 M or higher, i.e. 
much higher than the mass of the accreted shell M accr ~ 5.2 x 
10 7 M Q , inferred from UV and other observations during qui- 
escence. Therefore, T Pyx ejected far more material than it has 
accreted. 

There is no way to reconcile this discrepancy given the small 
uncertainty in the value of M; even if allowance is made for an 
uncertainty of a factor two, one obtains an upper limit to M accr 
that is smaller by a factor of at least ten than the theoretical value 
of M e j. Only for accretion rates higher than 4 xl0~ 7 M yr -1 
would the accreted mass M flccr be comparable with the estimated 
ejected mass M e j. However, these high rates would correspond 
to the steady-burning regime, while our detailed discussion of 
Sect. 14 definitely excluded this possibility. Further confirmation 
of our considerations can be found in the the very recent results 
of XMM observations that exclude the presence of a super-soft- 
source in T Pyx. 

The important point is that far more material appears to have 
been ejected during the last outburst of T Pyx than has been ac- 
creted by the white dwarf. This raises several doubts about the 
common assumption that the white dwarf in recurrent novae in- 
creases in mass toward the Chandrasekhar limit, and about the 
possible role of RNe as progenitors of SNIa. We note that Delia 
Valle and Livio (1996), based on statistical considerations on the 
frequency of occurrence of RNe in M3 1 and LMC, deduced that 
RNe are not a major class of progenitors of Type la supernovae. 
The behavior of T Pyx represents observational confirmation of 
this conclusion. 

Further confirmation for other RNe is required, and this high- 
lights the crucial need for accurate determinations of the most 



critical parameters of RNe, i.e. the mass accretion rate, and the 
mass and chemical composition of the shell ejected in a single 
outburst. 

In the case of T Pyx, at present, useful information can be ob- 
tained from highly spatially resolved spectrophotometry of the 
nebula that resolves its innermost part (associated with the last 
eruption), whose apparent radius should be by now larger than 
1". At the same time, spatially resolved observations of the outer 
portions of the nebula will shed light on its poorly known chem- 
ical composition and on its complex velocity structure. 
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